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SECTION  1 


UV  SOURCE  CHARACTERISTICS 
(Modification  of  Section  4  of  Volume  17) 


Subroutine  EUXFIT  did  not  require  a  change  for  ROSCOE-IR  per 
se,  but  we  have  done  so  anyway  so  that  the  routine  would  conform  to 
Sappenfield1 s  1975  report  [Sa-75a]  instead  of  the  1971  report  [FS-71a] 
by  Fajen  and  Sappenfield.  The  upgrading  involves  a  modified  formula 
as  well  as  altered  coefficients.  We  also  correct  an  error  near  the  end 
of  the  routine  where  the  total  hydro  yield  was  used  instead  of  that 
part  of  it  left  after  the  loss-cone  and  ion-leak  portions  had  been 
removed . 

Equation  (26)  in  Volume  17  is  still  used  with  two  exceptions: 

1.  For  the  second  energy  group,  use  x  =  1/K  instead 
of  x  =  In  K. 

2.  For  K  <40,  the  coefficients  in  Equation  (26)  for  the 
first  and  third  energy  groups  are  given  by  a 
different  set  of  values. 

As  implied  in  Volume  17,  we  again  add  the  coefficients  for  the  fifth 
and  sixth  energy  groups  given  in  Sa-75a  to  get  the  coefficients  we  use 
for  the  fifth  energy  group  in  the  code. 
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SECTION  2 

PARTITIONING  OF  THE  DEPOSITED  DEBRIS  ENERGY 
AND  CONSEQUENT  SPECIES  (SUBROUTINE  HPCHEM) 

(Revision  of  Section  7  of  Volume  17) 


2-1  INTRODUCTION 


The  principal  motivation  in  revising  Subroutine  HPCHEM  for 
ROSCOE-IR  has  been  to  provide  for  the  additional  species  carried  in  the 
high-altitude  grid  chemistry.  The  revision  is  extensive  yet  minimal. 

We  retain  the  original  philosophy  of  the  model  (based  on  a  light-energy 
deposition)  and  employ  a  reasonable  procedure  to  avoid  malfunction  when 
the  model  is  applied  to  heavy  energy-deposition  cases.  We  do  this  by  a 
modified  procedure  for  limiting  the  total  energy  deposited  so  as  to 
avoid  complete  burnout  of  diatomic  molecules.  We  also  correct  a  state¬ 
ment  in  Volume  17  for  the  maximum  allowed  elas tic-energy  deposition 
(HPEMAX) . 


Subroutine  HPCHEM  partitions  the  total  energy  lost  in  each 
cell  by  loss-cone,  ion-leak,  and  CHEX  particles  undergoing  inelastic 
and  elastic  collisions.  The  energy  is  partitioned  into  heavy-particle 
thermal,  electron  thermal,  dissociation,  ionization,  excitation,  and 
radiation.  A  model  of  heavy-particle  prompt  chemistry  is  also  used  to 


/  o  o 

compute  the  changes  in  number  densities  of  N(4S),  N(ZD),  N(ZP) ,  N 
N02,  0,  02(X3Z“),  02(a1Ag),  03,  CO,  C02,  H,  OH,  H02,  H20,  He, 


N0+,  0+,  ot.  and  electrons. 


N+,2N2 


NO, 


One  of  the  simplifications  made  in  depositing  the  heavy- 
particle  energies  is  that  the  range-energy  relation  used  in  Subroutines 
DEPO  and  CHXDEP  is  based  on  the  total  mass  density  of  a  cell  whereas  in 
Subroutine  HPCHEM  only  the  neutral  species  are  considered  as  targets. 
The  modeling  of  the  associated  heavy-particle  prompt  chemistry  is  based 
on  a  speculative  treatment  [HM-72b,  Section  6;  HL-73b,  Section  5], 
modified  by  the  need  to  interface  with  the  late- time  grid  chemistry. 
(The  rationale  for  the  heavy-particle  prompt -chemistry  modeling,  to  the 
extent  it  exists,  will  not  be  repeated  here.)  In  this  late- time  chem¬ 
istry,  for  atomic  oxygen  only  0(3D)  is  accounted  for  (implicitly)  as  an 


1 


V 


1 

A' 


r 

* 


excited  term.  Thus,  0(^S)  is  purposely  excluded  as  a  possible  reaction 
product  of  the  heavy-particle  prompt  chemistry. 

In  Tables  2-1  and  2-2  we  summarize  the  reactions  from  in¬ 
elastic  (ionizing)  and  elastic  (non- ionizing)  collisions  that  are 
modeled  in  the  prompt  chemistry  from  heavy-particle  energy  deposition. 
The  species  density  in  the  right-most  columns  are  initial  and  not  (more 
correctly)  current  values . 

2-2  ESTABLISHMENT  OF  PARAMETERS 

In  Subroutine  HPCHEM,  before  dealing  explicitly  with  the  par¬ 
titioning  of  the  species  and  energies  from  the  inelastic  and  elastic 
collisions,  we  must  take  a  number  of  preliminary  steps  and  establish 
several  parameters  used  in  the  modeling. 

1.  From  the  extended  core  storage  for  the  cell  of 
interest,  obtain  the  total  mass  density  p(g/cm^) 
and  the  number  densities  (cm” 3)  of  the  neutral 
species . 

2.  For  both  the  inelastic-  and  elastic-collision 
contributions,  sum  the  specific  energies  deposited 
by  the  loss-cone  (PEI  and  PEE) ,  ion-leak  (QEI  and 
QEE) ,  and  CHEX  (CEI  and  CEE)  particles. 

ei(erg/g)  =  PEI  +  QEI  +  CEI 
ee(erg/g)  =  PEE  +  QEE  +  CEE  . 

3.  Compute  the  number  density  of  nuclei  in  the  neutral 
species . 

A1  =  [N(4S)]  +  [N(2D>]  +  [N ( 2P) ]  +  [0]  +  [H]  +  [He] 

A2  =  [N2]  +  [NO]  +  [02(X)]  +  [02(a)]  +  [CO]  +  [OH] 

A3  =  [N02]  +  [03]  +  [C02]  +  [H02]  +  [H20] 

_  3 

n  =  neutral-particle  nuclei  cm 
=  A1  +  2  x  A2  +  3  x  A3  . 
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Table  2-1.  Heavy-particle  prompt -chemistry  reactions 
from  ionizing  collisions. 


Reaction 


N(4S)  •*  N+  +  e 
N(2D)  ->  N+  +  e 
N(2P)  -v  N+  +  e 


Branching  Number  of  Ionizing 


(n(4S) 


N*  +  e 
+  N+  +  e 
N0+  +  e 


NCS)  +  0+  +  e 
0(3P)  +  N+  +  e 
NO^  +  e  -*  NO  +  0 
NO  +  0+  +  e 
0(3P)  +  N0+  +  e 
0+  +  e 


( 0  2  +  ® 

I 0(3P) 


+  0+  +  e 


03  +  e  ->■  02  +  0 
02  +  0+  +  e 
0(3P)  +  02  +  e 


Ratio* 


Reactions,  cm" 


[NCS)][e]/n 

[N(2D)][e]/n 

[N(2P)][e]/n 

2[N2][e]/n 


2[N0][e]/n 


3[N02][e]/n 


[0][e]/n 


2[02][e]/n 


3[03] [e]/n 


a  =  1/1.35,  b  =  0.35a, 


(continued) 


n  =  number  density  of  nuclei  in  neutral  particles, 


Table  2-1  (Cont'd) 


Reaction 


Branching 

Ratio* 


(  C0+  +  e  C  +  0 

a 

CO 

, 

(  0(3P) +  C+  +  e 

b/2 

1  +  c 

C  +  0+  +  e 

b/2 

| 

(  C02  +  e  CO  +  0 

a 

co2 

-*•  < 

1 

[  0(3P)  +  C0+  +  e 

1  ~C+~ 

b/2 

CO  +  0+  +  e 

b/2 

H 

-4 

H+  +  e  -*■  H 

/ 

1.0 

0H+  +  e  -*•  0  +  H 

a 

OH 

'  0(3P)  +  H+  +  e 

1 

b/2 

1 

+ 

^  H  +  0  +  e 

b/2 

H02  +  e  -«•  OH  +  0 

a 

0(3P)  +  0H+  +  e 

0+H 

b/4 

HO  2 

OH  +  0+  +  e 

b/4 

H  +  02  +  e 

b/4 

02  +  H+  +  e 
<  H 

f 

b/4 

( 

H20+  +  e  ->-0H  +  H 

a 

h2o 

4  t 

1  H  +  0H+  +  e 

1  0+H 

b/2 

v  OH  +  H+  +  e 

H 

b/2 

He 

4- 

He+  +  e  ->  He 

1.0 

Number  of  Ionizing 
Reactions,  cm- 3 

2[C0][e]/n 

3[C02][e]/n 

[H][e]/n 

2 [OH] [e] /n 


3[H02][e]/n 


3[H20][e]/n 


[He][e]/n 
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Table  2-2.  Heavy-particle  prompt-chemistry  reactions  from 
non-ionizing  collisions. 


Reaction 


Branching  Number  of  Fast 

Ratio*  Cascade  Particles,*  cm~3 


Slow 


f 

4 


N 


N(4S)j 

+  N(2D) 

I  ? 

+  N(  r) 


N. 


NO 


NO. 


N(4S)  + 


N(4S)  + 


0(JP)  + 


NO  + 


0(JP)  + 


Fast 

N(4S) 

N(2D) 

N(2P) 

N(4S) 

N(2D) 

N(2P) 

0(3P) 

O^D) 

N(4S) 

IM(2D) 

N(2P) 

0(3P) 

0(XD) 

NO 


4/20 

10/20 

6/20 

4/20 

10/20 

6/20 

c/2 

d/2 

(4/20) /2 
(10/ 20) /2 
(6/20) /2 
c/2 
d/2 
1/2 


1 


[N]nc/n 


2[N2]nc/n 


2[NO]nc/n 


3[N02]nc/n 


(continued) 


* 


c  =  9/14,  d  =  5/14, 

nc  =  number  density  of  particles  in  cascade  of  elastic-collision 
energy . 
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Table  2-2.  (Cont'd) 


Branching  Number  of  Fast 

Reaction  Ratio*  Cascade  Particles,*  cm 


Slow 

Fast 

i  0(3P) 

c  | 

0 

/ 

)  1 

[0]n  /n 

(  0(XD) 

d  1 

c 

(02(X)  ) 

j  0(3P) 

c  j 

02  = 

+  ocp) 

4- 

I  2[02]n  /n 

(+O2 (a)) 

(  0(XD) 

d  1 

Z  C 

| 

( 

)  0(3P) 

c/2 

j 

°3 

1 

-  < 

1  °2 

| 

+ 

(  O^D) 

-  ; 

/  3[03]nc/n 

1 

[  0(3P) 

+ 

°2 

1/2 

> 

) 

(  0(3P) 

c/2  1 

I 

1  C 

+ 

) 

f 

CO 

< 

f 

(  0(iD) 

d/2 

>  2[C0]nc/n 

i 

[  0(3P) 

+ 

C 

1/2  j 

1 

1 

r 

(  0(3p) 

c/2  1 

| 

1 

I  CO 

+ 

<  -I 

f 

co2 

+  < 

1 

j 

(  0(XD) 

d/2  ’ 

>  3[C02]nc/n 

1 

0(3P) 

+ 

CO 

1/2  | 

| 

H 

-> 

\ 

H 

1.0 

[H]nc/n 

r 

(  0(3P) 

c/2  "j 

I 

I 

H 

+ 

1 

f 

OH 

-y  j 

1 

(  0(XD) 

d/2  J 

>  2[0H]nc/n 

1 

1  0(3P) 

< 

+ 

H 

1/2  J 

1 

(continued) 

Table  2- 


(Cont ' d) 


Number  of 
Cascade  Particles,*  cm 


-3 


Reaction 


Branching 

Ratio* 


4. 


Compute  the  threshold-energy  parameter  used  in  modeling  the 
cascade  of  elastic  energy. 

[N]  =  [N(4S>  3  +  [N(2D)]  +  [N(2P)] 

[02]  =  [02(X  3z-)]  +  [02(a  ^g)] 

Et  =  { 28 . 44[N]  +  39.04[N2]  +  26.04[NO]  +  12.48[N02]  +  23.4[0] 

+  20.48[02]  +  4.20[03]  +  44.44[C0]  +  21.80[C02] 

+  22 . 0[K]  +  17 . 60[OH]  +  9.40[H02] 


+  20.48[H20]  +  40.00[He]}/n  . 

In  this  equation  the  coefficient  of  [0]  includes  the 
replacement  of  em(0)  by  e_.(0)  +  er(0)  in  Eq .  (68)  of 
HM-72b  (or  Eq.  (46)  of  HL-73b) •  A  similar  remark 
applies  to  the  other  atomic  species.  Values  of  the 
parameters  used  are  given  in  Table  2-3. 

Compute  the  weighted  value  of  the  energy  removed  per 
collision  in  the  cascade  of  elastic  energy. 

e (eV)  =  { 14 . 22[N]  +  47.96[N2]  +  29.94[NO]  +  21.64[N02] 
+  11 . 70[0]  +  31.64[02]  +  23.5[03]  +  45.62[C0] 

+  34.30[C02] +  11[H]  +  30 . 20[OH]  +  23.40[H02] 


+  30.24[H20]  +  20[He] }/n 


6.  It  is  necessary  to  limit  the  deposited  specific  energies  to 
values  less  than  those  required,  according  to  the  model,  for 
complete  ionization  and  (effective)  dissociation  of  the 
neutral-particle  densities.  In  establishing  such  a  limit 
one  needs  to  recognize  that  the  model  used  here  for  parti¬ 
tioning  the  deposited  debris  kinetic  energy  is  a  linearized 
theory  and  hence  more  appropriate  for  a  light  than  heavy 
deposition  of  energy.  The  probability  for  destruction  of 
an  ambient  target  particle  is  assumed  to  be  proportional  to 
the  initial  (and  not  a  current)  density  of  the  particular 
target  species.  This  probability  is  also  weighted  by  the 
number  of  nuclei  in  the  target  particle.  That  is,  the  cross 
section  for  destruction  is  assumed  to  be  proportional  to 
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Table  2-3.  Parameters  for  the  idealized  cascade  of 
elastic-collision  energy  and  ionization 
energy  (eV) . 


Particle 

£  j  £__ 

a  m 

er 

ed+em+er 

ei 

N 

- 

2 . 22e 

12 

14.22 

14.53 

n2 

9.76 

2 . 22e 

12 

23.98 

15.58 

NO 

6.51 

1.46g 

7 

14.97 

9.27 

NO  2 

3.12 

U-i 

0 

O 

O 

» 

O 

7 

10.82 

9.78 

0 

- 

U-l 

O 

r-. 

0 

11 

11.70 

13.62 

°2 

5.12 

0 . 70f 

10 

15.82 

12.06 

°3 

1.05 

0 

0 

O 

1 

CsJ 

O 

s-x 

10 

11.75 

12.80 

CO 

11.11 

0 

O 

11 

22.81 

14.01 

co2 

5.45 

✓“N 

O 

1 

O 

N-/ 

O 

O 

11 

17.15 

13.77 

H 

- 

- 

11 

11 

13.60 

OH 

4.40 

0 . 70f 

10 

15.1 

12.94 

(2.71 

(HO  -  0) | 

ho2 

{l.99 

(H  -  02)[  0 . 35h 

9 

11.7 

11.53 

2.35  ) 

h2o 

5.12 

(H  -  OH) 

10 

15.12 

12.62 

He 

- 

- 

10 

20 

24.59 

e  From 

(10/20)  2 

: .  37  +  (6/20)  3.56  = 

2.22. 

^  From 

(5/14)  1. 

96  =  0.70. 

8  From 

the  mean 

of  the  values  for  e 

and  f. 

From  1/2  of  the  value  for  f. 


u 

A 
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the  number  of  atoms  in  the  target  particle.  Hence,  diatomic 
and  triatomic  particles  are  destroyed  more  quickly  than  the 
atomic  particles.  The  combination  of  this  feature  and  the 
linearized  theory  would  allow  the  diatomics  and  triatomics 
to  be  over-burned  out  if  the  deposited  energy  were  limited 
on  an  energetic  basis  alone,  i.e.,  if  one  were  to  use  the 
limit 

[e]  +  nc  <  NHp 

where  [e]  is  the  electron  density  produced  by  the  inelastic 
collisions,  nc  is  the  number  density  of  particles  in  the 
cascade  of  elastic  energy,  and  Nhp  is  the  number  density  of 
ambient  neutral  particles.  With  this  limit  the  diatomics 
and  triatomics  would  be  over-burned  out  by  the  factors 
NRp/(n/2)  and  NRp/(n/3),  respectively.  We  have  elected  to 
limit  the  energy  deposition  to  avoid  burnout  of  the 
diatomics,  to  prevent  certain  difficulties  in  Subroutine 
PCHEM  which  occur  if  the  initial  value  of  [ N2 3  is  zero. 

By  barely  preventing  the  burnout  of  diatomics,  we  do  allow 
the  over-burnout  of  triatomics,  which  we  correct  by  setting 
the  triatomic  species  number  densities  to  zero.  However, 
this  error  is  negligible  since  the  triatomic  species  number 
densities  are  very  small  in  the  altitude  region  in  which 
these  effects  occur. 

The  other  possible  procedures  include  (1)  using  target- 
independent  cross-sections  (which  seems  unrealistic)  or 
(2)  avoiding  linearized  theory  by  integrating  the  appro¬ 
priate  differential  equations.  Since  the  model  was  already 
complete  except  for  the  feature  of  limiting  the  deposited 
energy,  neither  of  these  two  possibilities  was  attractive. 

We  now  summarize  our  procedure  for  limiting  the  energy 
deposition. 


Define 


NHP(cm  '  =  A1  +  A2  +  A3 


'iraax 


~eraax 


=  28.5  e-|NHp  erg/g 
=  |  (Et  +  e)e1NHp  erg/g 


where  the  conversion  factor  e ^ 


- 1  3-1 

e^erg  eV  cm  g  )  =  eQ/p 


eQ  =  1.60  x  10"^  erg/eV  . 


1  O 

lo 


The  total  mass  density,  p,  is  appropriate  here  because  c. 
and  are  computed  for  the  total  mass  density.  The  1 
fractions  of  Njjp  destroyed  by  inelastic  and  by  elastic 
processes  (on  an  energetic  basis,  alone)  are,  respectively, 

®i  -  ei^eimax 


and 


I 


®e  ^e^emax 


If 


2-3 


F  =  NHp(g.  +  ge)/(n/2)  >  0.99  , 
one  must  reduce  and  by  the  factor  0.99/F. 
INELASTIC  COLLISIONS 


We  are  now  ready  to  deal  explicitly  with  the  inelastic 
collisions . 

1.  Compute  the  electron  density: 

[e]  «  e±/ (28.5  . 

2.  Partition  the  charged  species  from  inelastic  collisons. 

[N+]t  =  {b-1[N]  +  2[N2]  +  [NO]}  b[e]/n 

[N^] i  =  2[N2]  a[e]/n 
[NO+]i  =  { 2 [ NO ]  +  1.5(b/a)[N02]}  a[e]/n 
[0+].  =  {[NO]  +  [CO]  +  [OH]  +  1.5([N02]  +  [03]  +  [C02]) 
+  b-1[0]  +  2[02]  +  0.75[H02]}  b[e]/n 
[0+]i  =  (2[02]  +  1.5(b/a)([03]  +  0.5[H02])}  a[e]/n  . 


3. 


3.1  Record  the  corresponding  neutral-particle  densities 
produced  by  dissociative  ionization  in  inelastic 
collisions : 


[NfS)]di  -  (2[N2]  +  [NO]}  b[e]/n 


A 
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[N0]di  =  [N02J1.5  b[e]/n 

[0]d.  =  {[NO]  +  [CO]  +  [OH]  +  1.5([N02]  +  [O3]  +  [C02]) 
+  2[02]  +  0.75[H02]}  b[e]/n 
[°2]di  =  {t°3^  +  0.5[HO2]}  1.5  b[e]/n 
[CO]d.  =  [C02]  1.5  b[e]/n 

[H]di  =  {[OH]  +  0.75[HO2]  +  1.5[H20]}  b[e]/n 
[OH]di  =  {0.5[H02]  +  [H20]}  1.5  b[e]/n 

[C]di  =  [CO]  b[e]/n  . 

3.2  Prescribe  prompt  recombination  for  those  ions  not 
carried  in  the  late- time  grid  chemistry,  as  shown 
in  Table  2-4. 

3.3  Record  the  neutral-particle  densities  corresponding 
to  these  recombinations. 

[N0]r  =  [N02]  3a[e]/n 

[0]r  =  {3([N02]  +  [03]  +  (l+b/4a)[H02]  +  (b/2a)[H20] 

+  (l+b/2a)[C02J)  + 2([0H]  +  [CO])}  a[e]/n 
[02]r  -  [O3]  3a[e]/n 
[C0]r  =  [C02]  3a[e]/n 

[H]r  =  { (2+b/a) [OH]  +  a-1(1.5  b[H02] + [H]) 

+  3(l+b/a)[H20]}  a[e]/n 
[0H]r  =  {[H02]  +  [H20]}  3a[e]/n 
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Table  2-4.  Prescribed  prompt  recombination  of  (minor) 
charged  species  not  carried  in  late- time 
grid  chemistry. 


Reaction 

Branching 

Ratio* 

Number  of  Reactions 

cm"  3 

NO^  +  e  -  NO  +  0 

a 

3[N02] 

[e]/n 

0^  +  e  O2  +  0 

a 

3[03] 

[e]/n 

j 

(  a 

2[0H] 

[e]/n 

0H+  +  e  -*■  0  +  H 

b/4 

3[H02] 

[e]/n 

| 

b/2 

3[H20] 

[e]/n 

j 

f  b/2 

1 

2  [OH] 

[e]/n 

4- 

|  b/4 

3[H02] 

[e]/n 

H  +  e  ->•  H  < 

|  b/2 

3[H20] 

[e]/n 

1 

l  1-0 

[H]  [e 

]/n 

KOJ  +  e  -*■  OH  +  0 

a 

3[H02] 

[e]/n 

H20+  +  e  h.  OH  +  H 

a 

3[H20] 

[e]/n 

C02  +  e  CO  +  0 

a 

3[C02] 

[e]/n 

4 

(  b/2 

3[C0?] 

[e]/n 

CO  +  e  C  4-  0 

1 

2[C0] 

[e]/n 

C+  +  e  +  C 

b/2 

2  [CO] 

[e]/n 

He+  +  e  -*•  He 

1.0 

[He]  [ 

e]/n 

*a  =  1/1.35,  b  =  0.35a. 
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[He]r  =  [He]  [e]/n 

[C]r  =  { (2+b/ a) [CO]  +  (3b/2a)[C02]>  a[e]/n  . 

3.4  Record  the  decrease  in  electron  density  corresponding  to 
these  prescribed  recombinations . 

[e]r  =  ~{3([N02]  +  [03]  +  (l+b/2a) ([H02]  +  [C02]) 

+  (1+b/a) [H20] )  +  (2+b/a)([0H]  +  [CO]) 

+  a  ^([H]  +  [He])}  a[e]/n 

3.5  Record  the  densities  of  the  molecular  ions  in  Table  2-4. 
[N0j]r  =  [N02]  3a[e]/n 

[0^]^  =  [03]  3a[e]/n 

[0H+]r  =  { 2[0H]  +  (3b/2a)(0.5[H02]  +  [H20])}  a[e]/n 
[H02]r  =  [H02]  3a[e]/n 
[H20+]r  =  [H20]  3a[e]/n 
[CO^]^  =  [C02]  3a[e]/n 

[C0+]r  =  { 2[C0]  +  (3b/2a) [CO] }  a[e]/n  . 

3.6  Compute  and  collate  the  changes  in  neutral  species 
densities  for  inelastic  collisons. 

<S[N(4S)].  =  -[N(4S)]  Ill  +  [N(4S)]di 

6[N(2D) ] .  =  -[N(2D)]  Ill 
1  n 

6[N(2P)]  =  -[N(2P)]  Ill 

1  n 

6[N2].  =  -2[N2]  Ill 
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6 [NOJ  ± 
6[N02]i 
6[0], 
6[02]i 

*1 °3h 

6[CO]i 

6[C02]i 

6[H]. 

<S[OH]i 

6[H02]. 

6[H20]i 

6[He]. 

6j;c]. 


-2[N0] 

[e]  + 
n 

[“hi 

+ 

[NO] 

-3[N02] 

[e] 

n 

-CO]  til  +  [o]di  + 

[0] 

r 

-2[02] 

[e]  + 
n 

[°2hi 

+ 

[o2] 

-3[03] 

[e] 

n 

-2[C0] 

[e]  . 
n 

[“hi 

+ 

[CO] 

-3[C02] 

[e] 

n 

-[H]  i£l  +  [H]di  + 

[H] 

r 

-2 [OH] 

[e]  + 
n 

[0H]di 

+ 

[OH] 

-3[H02] 

[e] 

n 

-3[H20] 

[e] 

n 

0.0 

[c]di  + 

[C]r 

• 

3.7  The  changes  in  charged  species,  so  far, 


6[0+]  =  [0+]i 

6[<>J]  =  [0+]. 

6[e]  =  [e]  +  [e] 

r 

3.8  Compute  the  sums  of  changes  in  densities  of  nitrogen, 
oxygen,  carbon,  and  hydrogen  nuclei  and  of  charged 
species  as  measures  of  conservation  errors. 

Z6Nt  =  6[N(4S)]i  +  6[N(2D)]i  +  6 [N(2P) ] ± 

+  6 [NO] t  +  6[N02]i  +  6 [N+] ±  +  6[N0+]i 
+  2{6[N2]i  +  6[N2]i> 

E60i  =  6  [NO]  t  +  6  [0]  i  +  6[C0]i  +  6  [OH]  ±  +  6[H20]jL 
+  6[N0+]i  +  [0+]i 

+  2{6[N02]i  +  6[02]i  +  6[C02]i  +  6[H02]i 
+  6[0+].}  +  36[03]i 
I6Ci  =  6[C0]i  +  6 [C02]i  +  6[C]i 
E6H.  =  6 [H]  i  +  6 [OH]  .  +  5[H02]i  +  26[H20]i 
E6[Q+-Q“]  =  6[N+]  +  <S  [  N  J  ]  +  6[N0+]  +  6[0+] 

+  6[02]  -  <$[e] 

Partition  the  energy  from  inelastic  collisions  into  four 
modes  and  the  fraction  of  energy  in  each  mode. 

Electron  Thermal  Energy: 

E(eV/cm^)  =  1.5{[e]  4-  [e]r> 

fKE  =  C1  E/,£i 
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Dissociation  Energy: 

Di(eV/ cm3)  =  {9.76[N2]  +  6.51[NO]  +  5.12[02] 

+  11.11[C0]  +  4.40[OH]  +  1.5x(3.12[N02] 

+  1.05[03]  +  5.45[C02]  +  2.35[H02] 

+  5.12[H20])}  2b[e]/n 

f D  =  C1  Di/ci  ' 

Ionization  Energy: 

For  those  ions  that  are  not  carried  in  the  late- time 
grid  chemistry  and  therefore  undergo  prescribed,  prompt 
recombination,  the  effective  'ionization  energy'  correspond 
ing  to  the  molecular  ions  is  the  dissociation  energy  of  the 
neutral  particle;  for  the  atomic  ions  it  is  zero.  These 
facts  are  reflected  in  the  coefficients  of  the  following 
expression. 

I(eV/cm3)  =  14.53[N+]i  +  15.58[N2Ji  +  9.27[NO+]i 

+  13.62[0+]i  +  12.06[02]i  +  3.12[N02]r 


+  1 . 05[03]r  +  l.ll[CO+]r  +  5.45[C02]r 
+  4.40[OH+]r  +  2.71[H02]r  +  5.12[H20+Jr 


f-r  =  e-,  I/e. 


I1(erg/g)  =  exI 
Radiation  Energy: 


fR  =  1  -  <fKE  +  fD  +  fI> 


Ri (eV/ cm j  =  fRei/ei 
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ELASTIC  COLLISIONS 


1. 


We  now  deal  explicitly  with  the  elastic  collisions: 

Compute  the  total  number  density  of  particles  from  the 
cascade  of  elastic  energy: 


nc(cnf3)  = 


4fKEee 


where 

Ei-  -  £/3 

fKE  =  Et  +  e 

2.  Partition  the  neutral  species  from  elastic  collisions. 

2.1  The  fast-particle  densities  from  the  elastic-energy 


cascade  are: 

[N]^  = 

[N]nc 

/n 

[N]f2  = 

(2[N2 

3  + 

[NO]) 

nc/n 

[N]f  = 

[N]fi 

+ 

Mf2 

[NO] f  = 

[N02] 

1. 

5nc/n 

[o]n  = 

[0]  n 

c/n 

[°]f2  = 

{[NO] 

+ 

2[02]  + 

[CO]  +  [OH]}  nc/n 

[°]f3  = 

{[no2 

]  + 

[°3]  + 

[C02]  +  0.5[H02]}  1.5  nc/n 

f— n 

o 

1 _ 1 

ht> 

il 

[0]fi 

+ 

[°] f 2  + 

[°]f3 

II 

i — i 

(N 

O 

i _ i 

{[03] 

+ 

0.5[H02 

]}  1.5  nc/n 

[C0]f  = 

[  cc^] 

1. 

5  nc/n 
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[H]f  =  { [H]  +  [OH]  +  0.75[H02]  +  1.5[H20]}  nc/n 

[OH] f  =  {0.5[H02]  +  [H20]}  1.5  nc/n 

[He]f  =  [He]  nc/n 

[c]f  =  [CO]  nc/n  . 

2.2  The  slow-particle  densities  from  the  elastic-energy 
cascade  are: 

[N] s  =  [N] f 2 

[NO] g  =  [NO] f 

[O]  s  =  ^°^f2  +  E°]f3 

[02]s  =  { [03]  +  0.5[H02]}  1.5  nc/n 
[C0]s  =  [C02]  1.5  nc/n 
[H] g  =  {[OH]  +  0.75[H02]  +  1.5[H20]}  nc/n 
[OH] g  =  {0.5[H02]  +  [H20]}  1.5  nc/n 
[C]s  =  [CO]  nc/n  • 

2.3  Compute  and  collate  the  changes  in  species  densities 
for  elastic  collisions. 

<$[N(4S)]e  =  ~[N(4S)  ]  ^  +  ^j[N]f  +  [N]g 

5[N(2D)]e  =  -[N(2D)]  ^  ^  [N]f 

6[N(2P)]e  =  -[N(2P>]  ^  [N]f 

6[N2]e  =  -2[N2]  ^ 
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6[C0]e  =  -2[C0]  ^  4  [ CO] f  +  [CO]s 
6[C02]e  =  -3[C02]  ^ 

S[H]e  =  -[H]  ^  +  [H]f  4  [H]s 
6[OH]e  =  -2[OH]  ^  4  [OH] f  4  [OH]g 
[♦  6[HO2]0  =  -3[H02]  ^ 

l  n 

[  $[H20]e  =  -3[H20] 

<S[He]e  =  -[He]  ^  4  [He]f 
«[C]e  =  [C]f  4  [C]g  . 

,  r 

I  2.4  Compute  the  sums  of  changes  in  densities  of  nitrogen, 

\  oxygen,  carbon,  and  hydrogen  nuclei  as  measures  of 

conservation  errors. 

L6Ne  =  6[N(AS)]e  4  S[N(2D)]e  4  6[N(2P)]e 
4  [NO] e  4  6[N02]e  4  26[N2]e 


,'W* 


E60e  =  6[N0]e  +  6[0]e  +  6[C0]e  +  6[0H]e 
+  6[H20]e  +  36[03]e 

+  2{6[N02]e  +  6[02]e  +  6[C02]e  +  6[H02]e} 

L6Cq  =  <5[C0]e  +  <5[C02]e  +  <S[C]e 

Z6He  =  5[H]e  +  6[0H]e  +  6[H02]e  +  26[H20]e  . 

Partition  among  allowed  states  the  (total)  atomic-particle 
densities  from  the  elastic-energy  cascade. 

[N(4S)]  =  ^  [N]f  +  [N]s 

[N(2D)]  =  $  [N]f 

[N(2P)]  =  [N]f 

[0(3P)]  =  ^  [0]f  +  [o]s 

[0(3D)]  =  ^  [0]^ 

Partition  the  energy  from  elastic  collisions  into  four  modes 
Heavy-particle  thermal  energy: 


K  ( e  V  /  cnr  )  =  f,^/^ 
K-^erg/g)  =  tjK 
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Dissociation  energy: 


De(eV/cmJ)  =  {9.76[N2]  +  6.51[NO]  +  5.12[02] 


+  11.11[C0]  +  4. 40 [OH]  +  1.5(3.12[N02] 
+  1.05[03]  +  5.45[C02]  +  2.35[H02] 


+  5.12[H20])}  2nc/n 


Excitation  energy: 


Xe(eV/cm3)  =  (2.37  x  +  3.56  ^)[N]f 
+  1.96  x  [0] £ 


Xel(erg/g)  =  elXe 


Radiation  energy: 

Re(eV/cm3)  =  {0.5(12[N]  +  11[0]  +  11[H]  +  20[He]) 

+  12[N2]  +  7[N0]  +  10[02]  +  11[C0] 

+  10[0H]  +  1.5(7[N02]  +  10[03]  +  11[C02] 
+  9[H02]  +  10[H20])}  2nc/n  . 

2-5  COMBINED  RESULTS  FROM  INELASTIC  AND  ELASTIC  COLLISIONS 

We  can  now  combine  the  results  from  inelastic  and  elastic 
collisions . 

1.  Add  the  contributions  from  the  inelastic  and  elastic  colli¬ 
sions  to  obtain  the  total  values  of  dissociation  and  radia¬ 
tion  energies. 

D(eV/cm3)  =  D£  +  Dg 
D-^erg/g)  -  ejD 


R(eV/ cm3)  =  \  + 

R^Cerg/g)  =  ejR  . 

2.  Compute  the  total  fractional  partitions  of  energy  for  the 
six  modes. 

fK  =  £1  K/(ei  +  ee) 

fE  =  el  E/(£i  +  £e) 

fD  =  £1  D/<£i  +  ee> 

fX  el  Xe^ ^£i  +  £e^ 
e 

fI  =  £1  I/(£i  + 

fR  =  £X  R/(&i  +  eg)  . 

The  fractional  error  in  energy  conservation  is 

F  =  fK  +  fE  +  fD  +  fX  +  fI  +  fR  "  1  • 

e 

3.  Compute  the  changes  in  species  densities  for  both  inelastic 
and  elastic  collisions. 

$[N(4S)]  =  S[N(AS)].  +  6[N(4S)]e 

6[N(2D)]  =  6[N(2D)].  +  6[N(2D)]e 

6[N(2P)]  =  6[N(2P)]i  +  6[N(2P)]e 

6[N2]  =  6[N2].  +  6[N2]e 

6  [NO]  =  6[N0]i  +  <5[N0]e 

6[N02]  =  6[N02]i  +  6[N02]e 
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6[0] 

= 

fi[0]i  +  6[0]e 

6[02] 

= 

«S[02].  +  6[02]e 

«[02(X)] 

= 

S[02]  [02(X)]/[02] 

6 [ O2 (a) ] 

= 

6[o2]  u  -  [o2(x)]/[o2]> 

6[03] 

= 

6[03].  +  6[03]e 

6  [CO] 

(5  [CO]  t  +  6[C0]e 

6[C02] 

= 

6[C02]i  +  6[C02]e 

S[H] 

= 

6[H].  +  6[H]e 

6  [OH] 

= 

6[0H].  +  6[0H]e 

6[H02] 

= 

6[H02]i  +  6[H02]e 

<s[h2o] 

= 

6[H20],  +  6[H20]e 

6  [He] 

+  6[He]e 

«[C] 

6[C].  +  6[C]e  • 

4.  Compute  the  sums  of  changes  in  densities  of  nitrogen,  oxygen, 
carbon,  and  hydrogen  nuclei  as  measures  of  conservation 
errors . 

E6N  =  E6N.  +  E6N 
i  e 

E60  =  160.  +  E60 

l  e 

E6C  =  E6C.  +  E6C 
i  e 

E6H  =  E6Hi  +  E6He  • 
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5.  Compute  the  sum  of  electron  thermal  and  0(1D)  excitation 
energies  per  newly-formed  electron. 


6. 


E  (eV/electron)  =  {E  +  1 . 96(5/14) [0] }  6[e]  . 

r  e 

Store  changes  in  cell  quantities  in  temporary  scratch  area 
(BUF2) . 


BUF2  (1)  =  6[NCS)] 
BUF2 (2)  =  6[N(2D) ] 
BUF2  (3)  =  <5[N(2P)] 
BUF2 (4)  -  6[N2] 
BUF2 (5)  -  6 [NO] 
BUF2(6)  =  6[N02] 
BUF2(7)  =  6 [0] 

BUF2 (8)  =  6[02(X)] 
BUF2 (9)  =  6[02(a)] 
BUF2 (10)  =  fi[03] 
BUF2 (11)  =  6 [CO] 
BUF2(12)  =  6[C02] 
BUF2 (13)  =  6 [H] 


BUF2 (14)  =  6 [OH] 
BUF2(15)  =  6[H02] 
BUF2(16)  =  6[H20] 
BUF2 (17)  =  6 [He] 
BUF2 (18)  =  6[N+] 
BUF2(19)  =  6[N2] 
BUF2 (20)  =  6[N0+] 
BUF2 (21)  =  6[0+] 
BUF2 (22)  =  6[02] 
BUF2(23)  =  6[e] 


BUF2 (24)  =  E 


pe 


BUF2 (25)  =  Kx 
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SECTION  3 
UV  DEPOSITION 

(Modification  of  Section  9  of  Volume  17) 

3-1  INTRODUCTION 

The  total  computation  of  UV  deposition  requires  six  major 
subroutines,  as  noted  in  Section  9  of  Volume  17.  Here  we  give  a  very 
brief  overview  and  guide  to  the  changes  in  those  subroutines . 

The  modifications  to  Subroutine  EUXFIT  have  been  described 

in  Section  1.  Subroutine  PLINE  now  distinguishes  between  N(4S)  an-d 
2 

N(  D)  in  the  zeroing  of  line  integrals  performed  in  Subroutine  PINT. 
Subroutine  PINT  now  distinguishes  between  N(  b)  and  N(ZD)  in  forming 
the  line  integrals  of  the  species  densities.  Subroutine  BEDGE  now 
distinguishes  between  N(  b)  and  N(ZD)  in  determining  the  radius  to 
the  B-edge.  Subroutine  PHEAT  now  provides  for  all  the  additional 
species  in  ROSCOE-IR  and  for  N(  D)  as  an  attenuator  of  Group-U 
fluences.  The  extensive  changes  to  Subroutine  PCHEM  are  described  in 
Section  A. 


3-2  SUBROUTINE  PLINE  (Modification  of  Section  9.1.8.C  of 

Volume  17) 

A  new  integral  (SNDINT)  has  been  introduced  to  pertain  to  a 
2 

newly-carried  species, N(  D) .  This  new  integral  is  to  be  distinguished 
from  the  old  integral  (SNINT)  which  used  to  pertain  to  nitrogen  atoms 
(N)  without  distinction  as  to  the  ground  term,  but  which  now  pertains 
to  N(4S). 
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SUBROUTINE  PINT  (Modification  of  Section  9 . 2  of  Volume  17) 
Replace  the  integral  for  SNINT  by  the  two  expressions 


SNINT 


[N(4S)]dr,  cm"2 


[N(2D)]dr,  cm'2 


SNDINT 


Over  a  ray  segment  between  SR  and  SRN,  the  corresponding  contributions 
to  these  two  integrals  aie  DSN4SI  and  DSN2DI,  respectively.  SPECIE(I) 
now  denotes  the  concentration  of  N2,  O2,  N(4S),  0,  and  N(ZD),  respec¬ 
tively,  for  1=1,  2,  3,  4,  and  5. 

On  page  89  of  Volume  17,  DRHOI  is  misspelled  as  DHROI. 

3-4  SUBROUTINE  BEDGE 

3-4.1  Modification  of  Section  9.3.1  of  Volume  17 

The  old  line  integral  SNU  of  the  density  of  species  N  has 
been  replaced  by  the  two  line  integrals  SNU  and  SNDU  for  the  species 
N(^S)  and  N(^D) ,  respectively. 

3-4.2  Modification  of  Section  9.3.2  of  Volume  17 

The  index  I  in  SPECIE (I)  now  denotes  the  species  N2,  O2, 
N(^S) ,  0,  and  N(^D),  respectively,  for  1=1,  2,  3,  4,  and  5. 

3-4.3  Modification  of  Section  9.3.3  of  Volume  17 

In  Figure  20  on  page  94  of  Volume  17,  SPECIE (3)  is  replaced 

by  the  sum  (SPECIE (3)  +  SPECIE (5))  in  computing  EVEN,  RHOSUM,  and 
XMUBAR.  On  page  96  there  is  the  additional  line  integral, 

SNDU  =  SNDINT  +  SPECIE(5)  *  DEL  -  DSNDI , 

and  the  sum  over  J  is  extended  from  J  =  1,4  to  J  =  1,5  in  computing 
XMUBAR. 

3-5  SUBROUTINE  PHEAT  (Modification  of  Section  9.4  of  Volume  17) 

The  sum  over  J  is  extended  from  J  =  1,4  to  J  =  1,5  in 
computing  SUMIN(L) ,  SUMOUT(L) ,  PEFBAR(L) ,  FLUX(L),  and  XMU(L) .  See 
formulas  in  Section  9.4.1  and  in  9.4.3  of  Volume  17  (pages  98,  99). 


SECTION  4 


PROMPT  CHEMISTRY  ASSOCIATED  WITH  GROUP-X 
AND  GROUP-U  ENERGY  DEPOSITION  (SUBROUTINE  PCHEM) 

(Revision  of  Section  10  of  Volume  17) 


4-1  INTRODUCTION 

Subroutine  PCHEM  has  been  extensively  revised  in  detail  but 
not  in  overall  sturcture  in  going  from  ROSCOE-Radar  to  ROSCOE-IR.  New 
provisions  in  Subroutine  PCHEM  include: 

1.  The  additional  species  for  ROSCOE-IR. 

2.  An  increase  to  25  from  12  inputs  from  the  BUF2  array,  cor¬ 
responding  to  the  altered  output  from  Subroutine  HPCHEM. 

3.  A  revised  procedure  for  computing  the  ionization  and  species 
densities  resulting  from  the  Group-X  energy  deposition 
(computed  in  Subroutine  PHEAT) .  Included  are  dissociative 
ionization  of  N2  and  02»  species  and  excited  states  resulting 
from  non-ionizing  collisions,  and  an  improved  treatment  of 
CO2  destruction  with  protection  against  abortion,  when  [N2] 

is  very  small.  The  other  nine  minor  species  are  also  treated. 

4.  A  revised  procedure  for  computing  the  ionization  and  species 
densities  resulting  from  the  Group-U  energy  deposition.  The 
excited-state  species  N(^D)  is  now  included  in  the  line 
integrals  (performed  in  Subroutine  PLINE)  affecting  the 
Group-U  fluences.  Group-U  absorption  cross-sections  are 
provided  for  each  of  the  17  following  listed  neutral  species 
that  were  originally  designated  to  be  included  in  the 
ROSCOE-IR  grid  chemistry, 


I 

Species 

I 

Species 

I 

Species 

1 

n2 

7 

NO 

13 

HO  2 

2 

02(X) 

8 

NO  2 

14 

h2o 

3 

N(4S) 

9 

02(a) 

15 

co2 

4 

5 

0 

N(2D) 

10 

11 

03 

H 

16 

17 

CO 

He 

6 

N(2P) 

12 

OH 
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However,  owing  to  the  limited  types  of  ions  carried  in  the 
grid  chemistry  (N"^,  N0+,  O^,  N  ,  0+)  ,  we  prescribe  prompt 
dissociative  recombination  for  molecular  ions  and  radiative 
recombination  for  atomic  ions  not  carried  in  the  grid 
chemistry. 

5.  An  increase  to  35  from  18  outputs  in  the  BUF1  array. 

As  in  ROSCOE-Radar ,  Subroutine  PCHEM  first  combines  (a)  the 
changes  in  species  densities  due  to  debris  deposition,  computed  in 
Subroutine  HPCHEM,  with  (b)  the  previous  species  densities  of  a  cell 
and  then  computes  the  ionization  and  species  densities  from  the  Group- 
X  and  Group-U  energy  depositions. 

Subroutine  PCHEM  still  contains  numerous  compromises  and  ad¬ 
justments  made  necessary  by  using  energy  deposition  models  that  are 
(a)  instantaneous  and  (b)  mismatched  with  the  late -time  chemistry 
modules,  although  numerous  improvements  have  been  made. 

4-2  DETERMINATION  OF  INITIAL  SPECIES  FOR  COMPUTING  THE  PROMPT 

CHEMISTRY  ASSOCIATED  WITH  GROUP-X  AND  GROUP-U  ENERGY 
DEPOSITIONS 

Before  computing  the  prompt  chemistry  associated  with  the 
energy  depositions  of  Groups  X  and  U,  we  must  combine  the  changes  in 
species  densities  due  to  debris  deposition  with  the  previous  species 
densities . 

The  changes  in  the  species  densities  and  electronic  thermal 
and  excitation  energy  are  stored  in  the  BUF2  array  by  Subroutine 
HPCHEM.  The  number  of  these  variables,  defined  in  Table  4-1,  has  in¬ 
creased  to  25  for  ROSCOE-IR  from  12  for  ROSCOE-Radar.  These 
quantities  are  appropriately  combined  with  those  from  the  BUF  array 
to  get  the  starting  conditions  for  the  Group-X  and  Group-U  energy 
depositions.  We  check  charge  balance  (by  taking  the  difference  be¬ 
tween  positive  ions  and  electrons)  after  combining  the  HPCHEM-computed 
changes  in  species  densities  with  the  previous  densities  of  the  cells. 


Table  4-1.  Inputs  to  subroutine  PCHEK  from  BUF2 

array  for  those  cell  quantities  updated 
by  Subroutine  HPCHEM  as  a  result  of 
debris  deposition. 


I  BUF2 (I) 


1 

Change 

in 

[N(4S)] 

9 

cm"^ 

2 

Change 

in 

[N(2D)] 

9 

cm*3 

3 

Change 

in 

[N(2P)] 

9 

“  3 

cm 

4 

Change 

in 

[n2] 

9 

cm  3 

5 

Change 

in 

[NO] 

1 

cm*3 

6 

Change 

in 

[M02] 

9 

-3 

cm 

7 

Change 

in 

[0] 

9 

cm*3 

8 

Change 

in 

[02(X)] 

» 

cm  3 

9 

Change 

in 

[02(a)j 

9 

-3 

cm 

10 

Change 

in 

[o3] 

9 

cm*3 

11 

Change 

in 

[CO] 

9 

cm*3 

12 

Change 

in 

1 — 1 
CM 

o 

o 

1 — J 

9 

-  3 

cm 

13 

Change 

in 

[H] 

9 

cm  3 

14 

Change 

in 

[OH] 

9 

-  3 
cm 

15 

Change 

in 

[H02] 

9 

-  3 
cm 

16 

Change 

in 

[h20] 

9 

—  3 

cm 

17 

Change 

in 

[He] 

9 

cm*3 

18 

Change 

in 

[N+] 

9 

cm"3 

19 

Change 

in 

[N+] 

9 

cm-3 

20 

Change 

in 

[N0+] 

9 

cm"3 

21 

Change 

in 

[0+] 

9 

cm"3 

22 

Change 

in 

[°2] 

9 

cm"3 

23 

Change 

in 

[e] 

9 

-  3 

cm 

24  Change  in  electron  thermal  energy  and  0(D) 
excitation  energy  per  newly-formed  electron, 
eV/electron 

25  Change  in  heavy-particle  thermal  energy, 
erg/g 
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4-3 


IONIZATION  AND  SPECIES  DENSITIES  FROM  GROUP-X  DEPOSITION 


4-3.1  Introduction 

The  computation  of  the  species  densities  from  Group-X  deposi¬ 
tion  was  described  in  Volume  17  as  being  done  very  simply.  The  current 
computation  is  quite  complex,  resulting  partly  from  treating  additional 
species  required  for  IR  and  partly  from  attempting  to  make  the  compu¬ 
tation  more  realistic  for  both  radar  and  IR  purposes. 

To  facilitate  relating  this  text  to  the  coding,  we  describe 
the  major  steps  in  performing  the  overall  calculation.  Some  steps  are 
simple  enough  that  no  further  explanation  is  required.  In  other  cases 
the  steps  are  more  subtle  and  deserve  explanatory  information.  Some¬ 
times  we  will  provide  such  information  as  we  proceed;  at  other  times 
we  will  refer  to  Appendix  A  for  more  extensive  background  information. 

4-3.2  The  Major  Steps 

1.  Sum  the  N,  0,  C,  and  H  nuclei  densities  for  later  use  in 
mass-conservation  checks. 

2.  Preset  certain  variables  pertaining  to  Group-U  deposition. 
(While  this  step  has  nothing  to  do  with  Group-X  deposition 
per  se,  it  provides  for  the  possibility  that  a  Case-2  cell 
follows  either  a  Case-1  or  a  Case-3  cell.  For  a  Case-2  cell, 
no  computation  is  made  of  species  resulting  from  Group-X 
deposition . ) 

3.  3.1  Use  the  specific  energy  absorbed  from  Group-X 

deposition  (DAIXX(erg/g) )  computed  in  Subroutine 

PHEAT  and  the  mass  density  (RK0(g/cm3))  to  get 
the  corresponding  energy  density  (DELEX(eV/cm3) ) . 

DELEX  =  DAIXX  x  RH0/(1.6  x  10" 12)  .  (1) 

3.2  Determine  the  fraction  (FNEUT)  of  the  absorbed  energy 
to  be  assigned  to  deposition  in  the  neutral-particle 
gas.  (Here  we  consider  only  the  species  N,  N2,  NO, 

0,  Oo,  N+,  N^,  NO  ,  0+,  and  O^.)  The  prescription 
we  adopted  to  compute  FNEUT  is  described  in  Section 
A-2.2.2,  where  FNEUT  is  denoted  by  fn.  One  evaluates 
Equation  (A15)  with  use  of  Equations  (A16)  and  (A17) 
and  the  parameters  in  Table  A-4. 
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4. 


4.1 


Replace  the  value  used  in  Volume  17  for  the  energy 
per  ion  pair,  35  eV  per  ion  pair,  by  a  computed 
value  (EIP)  that  depends  on  the  mix  of  selected 
neutral  particles  (s  =  N,  No,  NO,  0,  and  O2) .  as 
described  in  Section  A-2.2.4.  Use  Equation  (A18) , 

EIP  E  E.  =  Ef  E.  (s) ,  (2) 

ip  g  s  ipw’  v  ' 

where  the  fractions  f  characterize  the  mix  for  the 
particles, 

f  =  [s]/l[s]  .  (3) 

s 

The  parameters  E^(s)  are  given  in  Table  A-4. 

4.2  Compute  the  ion-pair  density: 

AIPX(ip/cm3)  =  IP  =  DELEX  x  FNEUT/EIP  .  (4) 

5.  5.1  Divide  all  the  neutral  species  carried  in  the  HAG 

into  three  groups.  The  species  in  the  first 
group  are  chosen  because  they  are  the  principal 
neutral  species  and  because  we  have  an  improved 
technique  for  computing  the  energy  deposition. 

The  third  group  is  designated  to  facilitate  a 
special  treatment  for  CO2.  Group  2  contains  the 
remaining  minor  neutral  species.  The  groups  are: 

Group  1:  N,  N2>  0,  O2 

Group  2:  NO,  N02,  03,  CO,  H,  OH,  H02,  H20,  He 
Group  3 :  C02 

Sum  the  densities  [Jk]  eac^  °f  the  three  groups 
and  also  sum  the  group-sums: 

SUMGPk  =  I  [Jk]  k  =  1,2,3  (5a) 

Jk 

3 

SMNEUT  =  E  SUMGPi.  .  (5b) 

k=l  k 

5.2  Set  the  maximum  allowed  ion-pair  density  in  the 
neutral-particle  gas,  AIPXMX,  by  considering  only 
the  Group-1  particles.  For  the  current  model,  and 
for  considerations  of  late-tirae  species,  we  need 
to  set  the  number  of  J]_-particle  destructions  (by 
collisional  ionization,  and  for  molecules,  also 
dissociative  ionization,  collisional  dissociation, 
and  collisional  excitation)  per  ion  pair  in  a 
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Jn -particle  gas  as  (J^)  (J^  =  N,  N2,  0,  O2) .  To  set 
these  quantities  we  need  to  introduce  the  branching 


these  quantities  we  need  to  introduce  the  branching 
ratios  (%+  and  Rn+,  as  in  Section  A-2.1)  for  dissocia¬ 
tive  ionization  of  N2  and  O2  by  electron  impact  and  the 
ratios  (Fd(No)  and  Fd(02)  >  as  in  Section  A-2.3)  of 
collisional  dissociation  to  collisional  ionization. 

The  quantities  D.  (J,),  used  in  estimating  AIPXMX,  can 
be  set  as  " 

D.p(N)  =  1  Dip(N2)  =  1  +  (l-RN+)Fd(N2)  (6a, b) 

Dip(°)  =  1  Dip(02)  =  1  +  (l-RQ+)Fd(02) •  (6c, d) 

Finally,  we  compute  AIPXMX  from  the  expression 


AIPXMX  =  0.9  E 
Ji 


CJlln  CJl3 


LJ1J0 

i  TO 


L°lJo 

SUMGP, 


where  0.9  is  an  ad  hoc  safety  factor  and  J-,  =  N, 

N2 ,  0 ,  and  02 . 

5.3  Initially  distribute  the  Group-X  ion-pair  density 

among  the  three  groups  in  proportion  to  the  densities 
SUMGPk , 


AIPXk  = 


AIPX  x  SUMGP- 
E  SUMGP,. 


For  Group-1  and  -3  particles  we  will  later  redistri¬ 
bute  the  ion-pair  density  AIPX^  +  AIPX2.  For  Group-2 
particles,  we  simply  proceed  as  in  Volume  17  by 
distributing  the  ion-pair  density  AIPX2  according  to 
the  relative  abundance  of  each  of  the  Group-2 
particles . 

Determine  the  altered  species  densities  from  Group-X  deposi¬ 
tion  in  the  Group-2  particles  according  to  the  prompt- 
chemistry  reactions  assumed  in  Table  4-2,  where  [e]2  =  AIPX2 
and  n2  =  SUMGP2.  Note  that  collisional  ionization  and 
dissociative  ionization  of  the  Group-2  (and  Group-3)  par¬ 
ticles  lead  to  many  ions  not  carried  in  the  KAG  chemistry. 

To  obviate  this  difficulty  we  prescribe  immediate  dissicia- 
tive  recombination  for  such  molecular  ions  and  radiative 
recombination  for  such  atomic  ions. 


Table  4-2.  Prompt -chemistry  reactions  assumed  during 
x-ray  deposition  for  Group-2  (minor 
species)  particles. 


Branching  Number  of  Ionizing 
Reaction  Ratio*  Reactions,*  cm-3 


1 

f 

1 

N0+  +  e 

a  ] 

] 

\  N(4S) 

+  0+  +  e 

b/2 

[NO][e]2/n2 

NO  -*■  < 

4. 

1  0 

+  N  +  e 

b/2 

1 

[  N(4S) 

+  0 

c  > 

j 

f  no2  + 

e  -*•  NO  +  0 

3  ] 

1 

|  NO 

+  0+  +  e 

b/2  | 

>  [N02][e]2/n2 

N02  -  S 

r 

0 

+  N0+  +  e 

b/2 

k  NO 

+  0 

c  ) 

°3  +  ' 

2  -*•  02  +  0 

> 

a 

| 

°2 

+  0+  +  e 

b/2  1 

I 

03  i 

'  L. 

b/2  | 

/  [03][e]2/n2 

0 

+  O2  +  e 

<  °2 

+  0 

c  , 

' 

C0+  + 

e  -►  C  +  0 

\ 

a 

0  + 

C+  +  e 

b/2 

0 

0 

4 

_ 

1  C 

C 

+  0+  +  e 

b/2 

'  [CO][e]2/n2 

I 

1  c 

+  0 

c  ; 

(continued) 

a  =  0.75,  b  =  0.25,  c  =  0.8, 

[e]2  =  x-ray  deposition  ion-pair  density  assigned  to  Group  2, 
n2  =  density  of  group-2  particles. 


Table  4-2.  (Cont'd) 


Branching 

Reaction  Ratio 


H+  +  e 

!  H 

1.0 

0H+  + 

e  ->  0  +  H 

a 

0 

+  H+  +  e 

b/2 

H 

H 

4-  0+  4*  e 

b/2 

0 

+  H 

c 

ho2  + 

e  -*■  OH  +  0 

a 

0 

+  0H+  +  e 

b/4 

0+H 

OH 

+  0+  +  e 

b/4 

H 

+  02  +  e 

b/4 

°2 

+  H+  +  e 

b/4 

H 

OH 

+  0 

c/2 

H 

+  02 

c/2 

h2o+  + 

e  -*■  OH  +  H 

a 

K 

+  0H+  +  e 

b/2 

0+H 

OH 

+  H+  +  e 

b/2 

H 

OH 

+  H 

c 

He+  + 

e  •+  He 

1.0 

Number  of  Ionizing 
Reactions,  cm“3 


[H] [e]2/n2 


[OH] [e]2/n2 


[H20][e]2/n2 


[H20][e]2/n2 


[He] [e]2/n2 
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The  value  for  the  dissociative  ionization  branching  ratio 
(b)  is  a  nominal,  mean  value  suggested  by  averaging 
R^+ =  0.192  and  R a  =  0.302.  The  collisional  ionization 
branching  ratio  (a)  is  just  1-b.  The  collisional  dissocia¬ 
tion  branching  ratio  (c)  is  suggested  by  averaging 
[(l-RN+)Fd(N2)  =  1.00]  and  [ (l-R0+)Fd(02)  =  0.614]. 

Before  writing  the  changes  in  Group-2  particle  densities  we 
introduce  several  new  constants: 


Fl  =  a  +  0.5b 

(9a) 

F2  =  AIPX2/SUMGP2  5  [e]2/n2 

(9b) 

F3  =  i  +  C 

(9c) 

F4  "  F1  +  c 

(9d) 

F23  =  F2F3  ’ 

(9e) 

52[N(4S)]  =  F2(0.5b  +  c)[N0]q 

(10a) 

62[NO]  =  F2{F4[N02]q  -  F3[NO]q} 

(10b) 

62[N02]  =  -  F23[N02]o 

(10c) 

<52CO]  =  F2{F4([N02]o  +  [03]o  +  [C0]o  +  [0H]o 


+  (Fx  +  0.5c)[H02]q 

+  (0.5b  +  c) [NO]  +  0.5b[H20]Q}  (lOd) 

62[02]  =  F2{F4[03]o  +  (0.25b  +  0.5c)[H02]q}  (lOe) 

62[03]  =  -F23[03]o  (lOf) 

5 2 [CO]  =  -F23[CO]o  (lOg) 

62[H]  =  F2(F3[0H]o  +  (0.75b  +  0.5c)[H02]Q 

+  (a  +  1.5b  +  c)[H20]o)  (lOh) 


45 


F3[0H]o 


«2[oh]  -  f2{f4[h20]o  - 

+  (a  +  0.25b  +  0.5c)[H02]  }  (lOi) 

62Cho2:  =  -F23[H02]o  (lOj) 

62CH20D  =  -F23tH2°^o  (10k) 

52[He]  =  0  (101) 

S2[N+]  =  F2  0.5b[N0]o  (10m) 

<$2[N0+]  =  F2(a[NO]o  +  0.5b[N02]Q^  (lOn) 

<S2[°+]  “  f2  0-5b{[NO]Q  +  [N02]q  +  [03]q 

+  [CO]Q  +  [OH]q  +  0.5[H02]}  (IOo) 

62[0j]  =  F2  0.5b{[03]Q  +  0.5[H02]q}  (lOp) 

62[e]  =  F2([N0]o  +  b([N02]Q  +  [03]q 

+  0.5([C0]o  +  [OH]q  +  [H02]))}  (lOq) 

<S2C C]  =  F23tc°]0  '  (10r) 


One  can  show  analytically  that  the  N,  0,  C,  and  H  nuclei  are 
conserved  for  the  Group-2  reactions. 

Treat  the  Group-1  and  Group-3  particles  by  a  procedure  in 
which  some  of  the  steps  are  interleaved. 

7.1  Perform  the  first  nine  of  fifteen  steps  in  treating  the 
Group-1  particles.  The  background  information  for 
treating  the  Group- 1  particles  is  given  in  Appendix  A. 
The  information  in  Table  4-3  will  be  useful  in  con¬ 
sidering  the  first  seven  steps  for  the  Group-1 
particles . 


i  J 
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Table  4-3.  Ion-production  reactions  for  Group-1  particles 
(N2,  C>2,  N,  0)  during  x-ray  deposition. 


N 

0  •> 


Reaction 


•4  +  e 

N  +  N+  +  e 


O2  +  e 
0  +  0+  +  e 


N+  +  e 


0+  +  e 


Branching 

Ratio 


1  -  V  ) 
V  \ 


Number  of  Reactions* 
cm-  3 


(S"*" ;  e,  N2)ip  IP 


-  V  1 

Ro+  | 

I 

j  (S  ;  e,  02)ip  IP 

1.0 

(S+;  N)ip  IP 

1.0 

(S+;  0).p  IP 

The  molecular  reactions  include  the  modeling  of  just  the 
electron  slowing-down  collisions.  The  atomic  reactions 
include  the  modeling  of  both  the  photoabsorption  event 
and  the  electron  slowing-down  collisions.  This  table 
does  not  include  the  small  number  of  additional  N+  and  0+ 
ions,  included  in  the  model,  which  result  from  the  photo¬ 
absorption  events  in  N2  and  C^. 


Group-1  Step  1.  Determine  the  fractional  mix  of 
species . 


the  Group -1 


=  LiJ/2J[i]  (ID 

i 

i  =  N,  ^2*  0*  °2  * 


Group- 1  Step  2.  Determine  the  number  of  ions  from  species 
IT,  (S+;i),  per  typical  (1400-eV)  x-ray  photon  absorbed  in 
the  given  neutral-particle  mix  by  using  the  relations 


4.  « A  '  v- 

(S'; i)  =  a.f.1  +  b  f .  1 


(12) 


with  b'  given  by  Equation  (A13)  and  a.,  b.,  a.,  and  $. 
by  Table  A-l.  1  i  1  1 

Group- 1  Step  3.  Sum  the  ions  (S+;i)  from  each  of  the 
species  i: 


ST  =  £(S+;i)  •  (13) 

i 

Group -1  Step  4,  Normalize  the  quantities  (S+; i)  to  a  per- 
ion-pair  basis: 

(S+;i)ip  =  (S+;i)/ST  .  (14) 

Group- 1  Step  5 .  Decermine  the  ion  species  per  ion  pair  with 
account  of  the  dissociative-ionization  fractions, 


1^+  =  0.192  and  RQ+  =0.302  .  (15) 

Group- 1  Step  5.1.  For  the  initial  photoabsorption  events  in 
N2  and  O2,  let  and  Pq+  be  the  number  of  If1'  and  0+  ions 

per  ion  pair,  respectively: 


PN+  2fN2^T  (16a) 

PQ+  =  2f0^/Sx  .  (16b) 
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Group- 1  Step  5.2.  For  the  subsequent  electron  slowing-down, 
let 

4.  + 

(S  ;e>^2^ip  =  sum  ^  and  ^2  ^ons  Per  i°n  Pai-r 


=  (S+;N2).p  -  (17a) 

+  +  + 

(S  ;e ,02)ip  =  sum  of  0  and  02  ions  per  ion  pair 

-  <S+-°2>ip  -  V  '  (17b) 

-f  4.  4-  4 

Group- 1  Step  5.3.  The  number  of  N  ,  0  ,  ^  and  C>2  ions  per 
ion  pair  are: 

(N+).p  =  (S+;N).p  +  PN+  +  (S+;e)N2).pRN+  (18a) 

(0+)ip  =  (S+;0)ip  +  P0+  +  (S+:e’°2)ipR0+  (18b) 

(N+)ip  =  (S+;e,N2).p(l-RN+)  (18c) 

(0j)ip  =  (S+;e,02).p(l-R0+)  .  (18d) 

Group- 1  Step  6.  Determine  the  newly-formed  ion-species 
densities  corresponding  to  the  ion-pair  density: 

6l[N+]  =  (N+).pIP  (19a) 

51[0+]  =  (0+).pIP  (19b) 

«1[Nj]  =  (N2)ipIP  (19c) 

5l[oJ]  =  (02) ipIP  .  (19d) 


One  may  verify,  by  use  of  Equations  (18),  (17),  (14),  and 
(13),  that  the  sum  of  Equations  (19a)  through  (19d)  is  just 
IP. 

Group- 1  Step  7.  Collate  the  corresponding  changes  in  neutral- 
particle  densities  resulting  from  the  ionizing  events. 
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(20a) 


6U[N]  =  IP[(S+;e,N2).pRN+  -  (S+;N).p] 

6U[0]  =  IP[(S+;e,02)ipR0+  -  (S+;0).p]  (20b) 

6u[N2]  =  -  IP[(S+;N2).p  -  0.5Pn+]  (20c) 

6li[02]  =  "  Ipt(s+;°2)ip  “  0.5P0+]  .  (20d) 

To  demonstrate  that  nitrogen  and  oxygen  nuclei  are  conserved 
in  this  treatment  of  ion  production  in  x-ray  deposition,  we 
write 

6 (N  nuc.)  =  <5 (nitrogen  nuclei) 

=  <51[N+]  +  26x[N2]  +  6li[N]  +  26li[N2]  (21a) 

and  use  Equations  (19a),  (19c),  (20a),  (20c),  and  (17a): 

6(N  nuc.)  =  IP[(N+).p  +  2(N2)ip]  +  ^.[N]  +  261.[N2] 

=  IP[(S+;N).p  +  Pn+  +  (S+;e,N2).pRN+ 

+  (S+;e,N2)ip(l-RN+)2 

-(S+;N).p  +  (S+;e,N2).pRN+ 

+  PN+  -  (S+;N2).p2] 

=  IP[PN+  +  <S+;e,N2).p  -  (S+;N2).p]2 

=  IP[0 . 0 ] 2  .  (21b) 

A  corresponding  set  of  equations  and  conservation  obtain  for 
oxygen  nuclei, 

Group-1  Step  8,  Record  the  additional  changes  in  neutral- 
particle  densities  due  to  collisonal  dissociation: 
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(22a) 


col lisional  ionizations 


collisional  dissociations 

W°1  '  2V°2>  <l-V>  <S+,e,02)ipIP 
«lni[N2]  -  -  0.56lnt[N] 

«lni[02]  -  -  0.5«lnl[0]  . 


(22b) 

(22c) 

(22d) 


Group  1  Step  9.  Sum  the  changes  in  N,  0,  N2,  and  02  densi¬ 
ties  due  to  the  ionizing  events  and  collisional  dissociation: 


S-JN] 

-  +  6l„i™ 

(23a) 

6X[0] 

*  6li[«  +  W°3 

(23b) 

6x[N2] 

*  + 

(23c) 

<5-l[02] 

-  «uto2]  +  slni[°2]  ■ 

(23d) 

7.2  Treat  the  Group-3  particle  according  to  the  basic  procedure 
described  in  Section  A-3.  Equation  (A33)  relates  the 
decrease  in  CC>2  density  to  that  of  N2-  The  exponent  in 
Equation  (A33) ,  F(a),  given  by  Equation  (A29) ,  depends  on 
the  mean  value  of  the  fractional  ionization,  estimated 


from  the  following  equations: 

ao  = 

[e]0/{[e]Q  +  SMNEUT } 

(24a) 

te]f  = 

[e]Q  +  IP 

(24b) 

af  = 

[e]f/{[e]f  +  SMNEUT  -  IP} 

(24c) 

a  = 

0 . 5 (aQ  +  af )  . 

(24d) 

The  total  ionization  allocated  to  Group-1  and  Group-3  is 
iteratively  distributed  between  Group-1  and  Group-3. 
Temporarily  save  the  initial  value  of  6^[N2],  ( 6 ^ [N- 3 ) 0 • 
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Evaluate  Equation  (A33) : 


<53[C02]  =  -[C02]o  [1  -  {([N2]o 

+  61[N2])/[N2]o}F(a)]  (25a) 

<5 3 [ c°2 Jm  =  -  63[C02]  .  (25b) 

Assume  that  a  fraction  l/(l+c)  of  the  decrease  in  CO ^ 
density,  is  due  to  ionizing  reactions. 

I3  =  63Cc02]m/(1+c)  •  (26) 

The  ion-pair  density  I3  exceeds  the  ion-pair  density 
originally  allocated  to  Group-3,  AIPX3 .  Assume  this 
excess  in  ion-pair  density,  AI,  can  be  made  up  at  the 
expense  of  the  much  larger  density  of  Group-1  ions: 


AI  =  I3  -  AIPX3  .  (27) 

The  new  value  of  AIPXi  is  the  old  value  multiplied  by 
the  fraction  F  given  by 


F  -  1  -  AI/AIPXX 


(28) 


The  iterated  value  of  the  decrease  in  N2  density  is 


6{[N2]  =  F(61[N2])q 


(29) 


This  procedure  can  be  repeated  several  times  to  obtain 
a  satisfactory  convergence . 

When  a  satisfactory  value  of  F  is  adopted,  one  imposes 
charge  conservation  to  compute  the  following: 


AI  =  AIPX1(1-F)  (30a) 

I3  =  AIPX3  +  AI  (30b) 


<5  3  [  c°2  3m  =  i3(1+c) 


(30c) 


7.3  Record  the  iterated  values  of  the  changes  in  Group-1 
species  densities  by  multiplying  each  of  the  original 
values  by  F. 

7.4  Complete  the  following  six  steps  in  treating  the 
Group-1  particles. 

Group-1  Step  10.  Partition  the  total  change  in  02 
among  the  X-  and  a-states  of  0^  by  scaling  in 
proportion  to  the  initial  densities. 


ox[02(X)]  =  61[02]  [02(X)]o/[02]o  (31a) 

61[02(a)]  =  {1  -  [02(X)]o/[02]o}  .  (31b) 

Group- 1  Step  11.  Record  the  new  densities  of  X-  and 
a-states  of  02> 

[02(X)]1  =  [02(X)]q  +  51[02(X)]  (32a) 

[02(a)]1  =  [02(a)]Q  +  61[02(a)]  . 


Group-1  Step  12.  Compute  the  collisional  excitation  of 
02 (X)  to  0^ (a) .  By  using  Equation  (A26)  and  (A27)  we 
have : 

x  [0]o/[02]o  (33a) 


y 


(02(a)/0p 


ion  pair 


(33b) 


|  4.8  ,  x  <  0.6328 

(3.486  x"3/4  ,  x  >  0.6528 


(33c) 


[02(a)]ce  =  y6l[02]  • 


(33d) 


Group- 1  Step  13.  Record  the  new  (tentative)  densities 
of  tne  X-  and  a-states  of  C>2, 

[02(X)]2  =  [02(X)]1  -  [02(a)]ce  (34a) 

[02  (a)  ]  2  =  [02  (a)  ]  ^  4-  [*®2^a^ce  * 


(34b) 


r 
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Group- 1  Step  14.  Limit  the  number  of  02(a)  excitations 
if  the  02(X)  density  would  be  depleted,  i.e.,  if 

A[02(X)]  =  0.9[02(X)]1  -  [02(a)]ce  <  0,  (35a) 


use 


[02(a)]ce  =  0 • 9[°2(X) ]f  (35b) 
to  recompute  [02 (X) 3 2  and  [02(a)]2  . 

Group- 1  Step  15.  Partition  the  N  atoms  into  the  ground 
and  two  excited  states  (actually,  terms).  (See  Section 
A-2.3.) 


<51[N(4S)]  =  0.471 

(36a) 

61[N(2D)]  =  0.353 

^[N] 

(36b) 

61[N(2P)]  =  0.176 

(36c) 

7.5  Record  the  densities  of  the  Group-3  product  species 

(see  Table  4-4)  associated  with  (a)  the  ionizing  events 
and  (b)  the  collisional  dissociation  events  and  then 
combine  the  densities  for  0  and  CO. 


63[0].  =  I3  (37a) 

63[CO]i  =  I3F1  (37b) 

63[0+]i  =  I3  0.5b  (37c) 

63[e]i  =  63[0+].  (37d) 

63[C].  =  A3rO+J.  (37e) 

63[CO]ni  =  I3C  (38a) 

63[0]ni  =  I3C 


(38b) 


Table  4-4.  Prompt -chemistry  reactions  assumed  during 
x-ray  deposition  for  the  Group-3  particle. 


Reaction 


Branching  Number  of  Ionizing 
Ratio*  Reactions,  cm- 3 


CO, 


CoJ  +  e 

-y 

CO  +  0 

a 

0 

+ 

C0+  +  e 

b/2 

"TTT 

CO 

+ 

0+  +  e 

b/2 

CO 

+ 

0 

c 

I 


3 


See  Table  4-2. 


(39a) 


«3tO]  "  63[0]i  +  63^ni 
63CCO]  =  63[CO].  +  <53[CO]n.  .  (39b) 

8.  Compute  the  increased  atom  densities  and  decreased  electron 
density  from  imposing  dissociative  recombination  of  those 
molecular  ions  (N^j»,  N0+,  Cr£)  produced  by  Group-X  deposition 
and  carried  in  the  HAG  chemistry. 

+  + 

The  total  changes  in  densities  of  the  molecular  ions  N£,  NO  # 
and  0^  from  Groups  1  and  2  are 

{^[Np  +  0},  {0  +  62[NO+]},  and  { 6 x [ 0 J 3  +  62[02]}.  The 


assumed  dissociative  recombination  reactions  are: 

N2  +  e  -  N(4S)  +  N(2D)  (40a) 

oj  +  e  =  0  +  0  (40b) 

N0+  +  e  =  0.25N(4S)  +  0.75N(2D)  +  0  .  (40c) 

Thus  the  increased  atom  densities  from  dissociative 
recombination  are: 

AN2[N(4S) :  =  6x[N2]  (41a) 

An  [N(2D)]  =  6x[Np  (41b) 

Anq[N(4S)]  =  0.2562[NO+]  (41c) 

Anq[N(2D)]  =  0.7562[NO+]  (41d) 

A  no  t 0  ]  =  <5  2  [NO"1"]  (Ale) 

AQ  [0]  =  2{61[02]  +  6 2 [Op  }  .  (41f ) 

The  change  in  electron  density  is: 

A[e]  =  -C61[N2]  +  62[NO+]  +  <S1CoJ]  +  «2[0p]  }  .  (41g) 
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Collate  all  species  densities  after  Group-X  deposition  and 
aforesaid  dissociative  recombination  of  molecular  ions: 


[N(4S)]xf  =  [N(4S)]0  +  a-^N^S)]  +  6N2[N(4S)] 

+  5 2 L N ( )  ]  +  AN0[N(4D)]  (42a) 

[N(4s)]xf  =  [N(2D)]q  +  61[N(2D)]+  AN2[N(2D)] 

+  Ano[N(2D)]  (42b) 

[N(2p)]xf  =  £N(2p)]0  +  5x[N(2P)]  (42c) 

[N2]xf  -  [N2]0  +  61LN23  (42d) 

[N0]xf  =  [N0]0  +  a2[NO]  (42e) 

[N02]xf  =  [N02]q  +  62[N02]  (42f) 

[0JX£  =  E 0 3 q  "b  5^0]  +  4^q[0]  +  Aq^LO] 

+  a2[0]  +  63[0]  (42g) 

C  °2 (x) ]xf  =  [02(X>]2  +  62[02]  (42h) 

[°2(a)]xf  =  [02(a) ] 2  (42i) 

[02JX£  =  [03]0  +  62[03]  (42j ) 

[c°]xf  =  [C0]0  +  62[CO]  +  6  3 [ CO  3  (42k) 

[C02]xf  —  [C02]q  ^ 3 E CO 2 D  (421) 

[ H] xf  =  [H3  0  +  52[H]  (42m) 

[0H]xf  -  [OH] 0  +  62[OH] 


(42n) 


(42o) 


£H02^xf  - 

[H2°Lf  - 

[H02]q  +  62[H02] 

[H20]q  +  62[H20] 

(42o) 

(42p) 

[HeLf  - 

[He]0 

(42q) 

- 

[N+]q  +  61[N+]  +  62[N+] 

(42r) 

t^Xf  - 

[n2]0 

(42s) 

[N°+]xf  - 

[N0+]0 

(42t) 

[°+]xf  =  [0+]0  +  S1[0+]  +  62[0+]  +  63[0+]  (42u) 

[°2]xf  =  [°2^0  (42v) 

E e 1 xf  =  [e]Q  +  AIPX  +  A[e]  -  {AIPX2  -  62[e]} 

-  {I3  -  .  (42w) 


10.  Sum  the  final  densities  of  N,  0,  C,  and  H  nuclei,  followed  by 
finding  the  fractional  increases  in  the  final  values  relative 
to  the  initial  values,  as  measures  of  conservation  errors  in 
treating  the  Group-X  deposition.  One  can  show  analytically, 
as  well  as  numerically  in  the  code,  that  the  nuclei  are  con¬ 
served.  To  test  charge  balance,  divide  the  difference  of  the 
final  values  of  positive  charge  density  and  electron  density 
by  the  electron  density. 

4.4  IONIZATION,  SPECIES  DENSITIES,  AND  RESULTANT  PRESSURE  FROM 

GROUP-U  DEPOSITION 

In  computing  the  species  from  Group-U  deposition,  we  consider 
three  cases,  just  as  we  did  in  Subroutine  PHEAT: 

Case  1:  B-edge  radius  (SRR(2))  is  smaller  than  the 
target-cell  entry  radius  (SR) . 

Case  2:  B-edge  (SRR(2))  is  greater  than  the  target¬ 
cell  exit  radius  (SRN) . 

Case  3:  B-edge  radius  (SRR(2))  is  between  the  target¬ 
cell  entry  radius  (SR)  and  the  exit  radius  (SRN) . 


58 


The  total  attenuation  factors  (FSIGU(J))  are  computed  from 

2 

the  summed  products  of  Group-U  fluences  (F(L)  =  FLUX(L) ,  photons/cm  ) 

2 

and  cross  sections  (SIGU(L,J),  cm  ).  (Recall  that  the  Group-U  fluences 
computed  in  Subroutine  PHEAT  are  effective  fluences,  consistent  with 
the  absorbed  energy  computed  as  the  difference  of  the  (steradianal)  en¬ 
ergy  entering  and  leaving  the  cell.)  The  total  attenuation  factors  are 

5 


A  ( J)  =  FSIGU(J)  =  Ef(L)  X  SIGU(L,  J)  ,  J  =  1,17  . 

L=1 

(43) 

The  attenuation  factors  for  ionization  per  se  of  N2 
and  02(a)  are  obtained  by  subtracting  terms  for  dissociation: 

,  02(X), 

X^l)  =  A  (1)  -  SIGU  (1 , 1)F (1) 

(44a) 

A.(2)  =  A (2)  -  SIGU(1 , 2) F(l) 

(44b) 

Xt(9)  =  A (9)  -  SIGU(1 , 9)F (1)  . 

(44c) 

The  total  survival  probability  of  species  J  is 

P(J)  =  PROBSU(J)  =  e'X(J)  . 

(45) 

The  branching  ratios  a-^,  a2,  and  a^  for  tit,  0^,  and 
production  from  N2 ,  02(X),  and  02(a),  respectively,  are 

°2 

aL  =  A.(l)/A(l) 

(46a) 

a2  =  Ai(2)/A(2) 

(46b) 

a9  =  Ai(9)/A(9)  . 

(46c) 

SIGU(L,J)  is  the  absorption  cross-section  of  Species  J  for 
Subgroup-L  photons,  given  in  Tables  4-5a,  4-5b,  and  4-5c.  The  cur¬ 
rently  adopted  values  are  in  the  columns  headed  Vol.  17-1  (and 
entered  as  data  in  Block  Data  BLOCKH)  and  may  be  compared  (where  there 
are  corresponding  species)  with  the  values  used  for  ROSCOE-Radar  in 
the  columns  headed  Vol.  17.  The  reactions  considered  for  Group-U 
deposition,  including  instantaneous  dissociative  recombination  of  all 
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A  2 

Table  4-5b.  UV  photoabsorption  cross-sections  (N(s) ,  N(^D)  , 
N(2P) ,  0,  NO) . 


ROSCOE 

Group 

Number 

UV 

Bench¬ 

mark 

Group 

Number 

a(N(4S)) 

o (N(2D) ) 

a(N(2P)) 

HL- 7  3b 

Vol . 
17-1 

Vol. 

17- 

HL-73b 

Vol. 

17-1 

HL-73b 

Vol. 

17-1 

\ 

( 

1 

0 

| 

4.5  ) 

4.85  1 

1  < 

i 

f 

\ 

2 

0 

f 

•  0 

0 

4.1  j 

’  4.3 

7.6  ( 

'  6.2 

( 

[ 

3 

10.2 

) 

9.95  \ 

8.6  j 

2  < 

1 

F 

( 

4 

10.7 

l 

( 

10.4 

10 

9.3  j 

•  9.7 

10.9  i 

>  9.8 

1 

f 

5 

10.7 

> 

| 

9.4  ' 

| 

10.9  S 

| 

1 

| 

6 

10.9 

1 

[ 

9.4  1 

10.9  1 

[ 

3  1 

1 

7 

11.0 

| 

>  10.9 

10 

10.6  | 

/10.0 

10.8  | 

>10.9 

( 

[ 

8 

11.0 

) 

10.6  . 

) 

10.8  ; 

| 

4  j 

\ 

9 

11.3 

\ 

i 

1  11.1 

10 

10.3  ! 

1 

>  9.7 

10.46  j 

i 

-  9.8 

1 

\ 

10 

10.9 

s 

1 

9.1  j 

1 

9.1  j 

5 

1 

11 

5.9 

l 

5.9 

10 

2.6  | 

2.6 

4.0  \ 

4.0 

ROSCOE 

Group 

Number 

UV 

Bench¬ 

mark 

Group 

Number 

o  (0) 

a  (NO) 

Vol. 

HL-73b  17-1 

Vol  . 
17d 

HL-73b 

Vol. 

17-1 

*  ! 

1 

2 

0  !■  ° 

0  1 

0 

10.3  1 

28.9  1 

20 

19 


21 
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Table  4-5c.  UV  photoabsorption  cross-sections  (N02,  OoC^A  ), 
03>  H,  H20,  OH,  H02,  C02,  He).  8 


ROSCOE 

a 

(N02) 

o(02(1 

y> 

o(03) 

o  (H) 

Group 

Number 

Vol. 

17-1 

Ref. 

Vol. 

17-1 

Ref. 

Vol. 

17-1 

Ref . 

Vol. 

17-1 

Ref. 

1 

20 

SH-66a 

2.6 

(g) 

'x.lO 

SH-66a 

0 

2 

19 

(f) 

20 

(g) 

23 

SH-66a 

5.3 

PK-76c 

.  3 

21 

(f) 

22 

(g) 

34 

SH-66a 

4.0 

PK-76c 

4 

21 

(f) 

24 

(g) 

27 

SH-66a 

1.8 

(i) 

5 

23 

(f) 

16 

(g) 

9 

(h) 

%0 

(i) 

ROSCOE 

a(H20) 

a  (OH) 

a(H02) 

o(C02) 

IVv  V  ViJ 

Group 

Number 

Vol. 

17-1 

Ref. 

Vol. 

17-1 

Ref. 

Vol. 

17-1 

Ref. 

Vol. 

17-1 

Ref. 

1 

9 

SH-66a 

9 

(k) 

7 

(1) 

14 

2 

16.5 

SH-66a 

17 

(k) 

18 

(1) 

21 

;  3 

16 

SH-66a 

16 

(k) 

19 

(1) 

35.8 

t  4 

14 

(j) 

14 

(k) 

19 

(1) 

31.5 

5 

7 

(j) 

7 

(k) 

11 

(1) 

21.5 

Footnotes  for  Tables  4-5a,  4-5b,  and  4-5c: 


r 

* 


4 


f 

* 


L 


(a)  Sa-75a  quotes  bin  edges  at  6.2,  13.6,  15.6,  17,  30,  and  82  eV. 

(b)  HL- 7 3b ,  Section  2.3.2. 


(c)  Multiplets  included  in  the  photodeposition  code  FOTOKEM  [HL-73a]; 
multiplet  number  from  WS-66. 


(d) 


Ch) 


(1) 


(m) 

(n) 


(P) 


Volume  17  of  The  ROSCOE  Manual,  quoting  Table  A-l  of  FS-73  (in 
turn  from  DD-70b) . 


(e)  Variation  in  cross  sections  for  multiplet  lines  is  large. 

O 

(f)  No  data  are  available  at  wavelengths  below  vlOOO  A;  therefore,  we 
assumed  cr(N02)  ^  a (NO) . 

(g)  No  data  are  available  for  09(^A  );  therefore  we  assumed 

o(09(1a„))  %  a(09).  2  S 


2k  ug"  “  "^2' 

Based  on  extrapolation  of  data  in  SH-66a. 


(i)  Based  on  extrapolation  of  data  in  PK-76c. 


(j) 


In  SH-66a,  data  between  195  and  1085  A  of  J.  Romand  (private 
communication)  are  given.  The  absolute  values  of  these  data  are 
significantly  higher  than  those  obtained  by  Watanabe  and  Jursa 
[1964]  or  Metzger  and  Cook  [1964a]  as  reported  in  Hu- 71b.  To 
obtain  a(H20)  listed,  the  data  of  J.  Romand  were  averaged  and 
scaled  to  the  values  of  Metzger  and  Cook. 

(k)  No  data  are  available  for  OH;  therefore  we  assumed 
a (OH)  ^  a(H20). 

No  data  are  available  for  H0-> ;  therefore  we  assumed 


aata  are  avarraDie  ror  tiuo ;  cnererore  we  assum< 
o(H02)  %  0 . 5[o(02)  +  a (OH) ]  $  0.5(a(02)  +  a(H20)] 

Derived  from  measurements  reported  in  SH-66a  and  LC-73. 

0  to  as  large  as  160  Mb  in  the 


Cross  sections  range  from  about 
wavelength  range  of  Group  1. 

(o)  Cross  sections  range  from  about  15  to  120  Mb  in  the  wavelength 
range  of  Group  2 . 


Cross  sections  range  from  about  15  to  60  Mb  in  the  wavelength 
range  of  Group  3. 
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molecular  ions,  are  given  in  Table  4-6. 

To  facilitate  writing  the  differential  equations  for  the 
changes  in  species  densities  during  the  Group-U  deposition,  we  intro¬ 
duce  the  following  notation. 


S(J)  =  A(J)/F  (47a) 

5 

=  E  SIGU(L, J)f (L)  (47b) 

L=1 

where  f(L)  =  fractional  spectral  fluence 
=  dF  (L)  /  dF 

dF(L)  =  differential  fluence  of  Subgr'oup-L  photons 

dF  =  differential  total  fluence. 

The  differential  equations  for  the  species  densities  during 
the  Group-U  deposition  are: 

d[N2]/dF  =  - [N2] S (1)  (48a) 

d[02(X)]/dF  =  -[02(X)]S(2)  +  [03]S(10)  (48b) 

d[N(2S)]/dF  =  -[N(4S)]S(3)  +  2g1[N2]S(l)  +  32[NO]S(7)  (48c) 

d[0] / dF  =  - [0] S (4)  +  2[02(X)]S(2)  +  [N0]S(7) 

+  [N02]S(8)  +  2[02(a)]S(9)  +  [03]S(10) 

+  [0H]S(12)  +  [H02]S(13)  +  [C02]S(15) 

+  [CO] S ( 16)  (48d) 

d[N(2D) ] / dF  =  ~[N(2D) ]S(5)  +  2(1  -  81)[N2]S(1) 

+  (1  -  02) [NO ]S(7)  (48e) 
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Table  4-6.  Reactions  for  Group-U  deposition. 


Species  UV 

Number  Reaction*  Subgroup-L 


2g]N(4S) 

+ 

2(1 

- 

&1)n(2d) 

1 

1 

n2 

->  j 

1  N^  +  e 

-> 

23]N(4S) 

+ 

2(1 

- 

61)N(2D) 

2.5 

1 

| 

0 

+ 

0 

1 

2 

09(X) 

-> 

1 

l  4- 

|0j  +  e 

-> 

0 

+ 

0 

2,5 

3 

N(4S) 

-> 

N+ 

+ 

e 

2,5 

4 

0 

-► 

0+ 

+ 

e 

2,5 

5 

N(2D) 

->■ 

N+ 

+ 

e 

1.5 

6 

N(2P) 

-> 

N+ 

+ 

e 

1,5 

7 

NO 

N0++  e 

-> 

62N(4S) 

+ 

(1  - 

-e2 

)N(2D)  +  0 

1,5 

8 

CN 

O 

53 

-> 

NO^  +  e 

-> 

NO 

+ 

0 

1,5 

1 

( 

0 

+ 

0 

1 

9 

07(a) 

1 

z. 

1 

(  Oj  +  e 

-+ 

0 

+ 

0 

2,5 

1 

( 

07(X) 

+ 

0 

1 

10 

°3 

J 

l°3  +  e 

-+ 

02(X) 

+ 

0 

2,5 

11 

H 

-V 

H+  +  e 

->■ 

H 

2,5 

( 

0 

+ 

H 

1 

12 

OH 

-> 

(  0H+  +  e 

0 

+ 

H 

2,5 

( 

OH 

+ 

0 

1 

13 

H0o 

) 

i  _L 

z 

(  HOj  +  e 

-+ 

OH 

+ 

0 

2,5 

(continued) 
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Table  4-6.  (Cont 'd) 


Species 

Number 


Reaction 


UV 

Subgroup -L 


■  ,,;£V  * 


d[N(2P)]/dF 

= 

-[N(2P)]S(6) 

(48f) 

d[NC]/dF 

= 

- [NO]S (7)  +  [N02]S(8) 

(48g) 

d[N02]/dF 

= 

-[N02]S(8) 

(48h) 

d[02(a)]/dF 

= 

-[02(a)]S(9) 

(48i) 

d[03]/dF 

= 

- [O3] S (10) 

(48j ) 

d[H]/dF 

= 

- [OH] S (12)  +  [H20]S(14) 

(48k) 

d[OH] / dF 

= 

-[0H]S(12)  +  [H02]S(13)  +  [H20]S(14) 

(481) 

d[H02]/dF 

= 

-[H02]S(13) 

(48m) 

d[H20]/dF 

= 

-[H20]S(14) 

(48n) 

d[C02]/dF 

= 

~[C02]S(15) 

(48o) 

d[CO]/dF 

= 

- [C0]S (16) 

(48p) 

d[He]/dF 

as 

0  . 

(48q) 

After  introducing  the  following  notation,  we  write  the  solu¬ 
tions  to  the  differential  equations. 

G(J,K)  =  G(X,P,J,K) 

=  PCFUN(A ,P, J,K) 

=  [P (K)  -  P(J)]X  (K)  / [  A  (J)  -  X  (K)  ]  (49) 
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H(I,J,K)  =  H(X  ,  P ,  I ,  J  ,  K) 


=  PCFUNC ( X , P , I , J , K) 


X (I)G(K , J)  -  X (J)G(K, I) 
X(I)  -  X(J) 


PCFUN  and  PCFUNC  appear  as  function  routines  in  the  code. 


[N2] 

[N9]  P(l) 

2  uf 

2  xf 

[02(X)] 

2  uf 

[09(X)]  P (2) 
2  xf 

+  [O3]xfG(2,10) 

[N(4S)]uv  = 

[N(4S)]xfP(3) 

+  23-,  [N9]  G(3,l)  +  B9[N09]  H(7 , 8, 3) 

1  2  xf  2  2  xf 


[°]uv  =  C0]xfp(4)  +  2[02(X>]xfG(4,2)  4-  [NO]xfG(4,7) 


+  [N09]  {G(4, 8)  +  H(7 , 8 , 4) }  +  2[02(a) ]  f G(4,9) 

xf 


+  [0,]  {G (4, 10)  +  2H(2 ,10,4)}  +  [0H]vrG(4, 12) 
J  xf  X 


+  [H09]  {G(4, 13)  +  H(12 , 13 , 4) }  +  [C0]vfG(4, 16) 

2  xf 


+  [C09]  {G(4,15)  +  H(16,15,4)}  +  [H90]  H(12,14, 

2  xf  xf 


[N(2D)]uv  =  [N(2D)]xfP(5)  +  (l-32)[NO]xfG(5,7) 

+  2(1-3,)[N9]  G (5 , 1)  +  (1- 89) [N09 ]  H(7 , 8 , 5) 

1  2  xf  2  2  xf 

[N(2P)]uf  =  [N(2P)xfP(6) 


(50) 

(51a) 

(51b) 

(51c) 


4) 

(51d) 

(51e) 

(51f) 
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(51g) 


[N0]uf  = 

[N0]vfP(7)  +  [N0?]  G(7 , 8) 

xt  z  xf 

(51g) 

[N0?] 

[N0?]  P (8) 

(51h) 

1  uf 

Z  xf 

[02 (a) ] 
z  uf 

[02(a)]xfP(9) 

(51i) 

[o3] 

[03]  P(10) 

(51j) 

J  uf 

J  xf 

For  hydrogen,  set  X  ( 1 1 )  =  0  and  P(ll)  =  1. 


[H]uf  =  [H]xf  +  [0H]xfG(11.12)  +  [H20]{G(11,14) 

+  H(14, 12, 11) }  +  [HO,]  H(13 , 12 , 11)  (51k) 

z  xf 


[°H]uf  = 

[OH]  ,P(12)  +  [HO,]  G(12 , 13)  +  [H,0]  G(12,14) 

x£  Z  xf  Z  xf 

(511) 

[HOp] 

[HOp]  P (13) 

(51m) 

1  uf 

Z  xf 

[HpO] 

[HpO]  P(14) 

(51n) 

1  uf 

Z  xf 

[co2] 

[COp]  P (15) 

(51o) 

Z  uf 

Z  xf 

[C0]uf  = 

[CO]  fP(16)  +  [COp]  G(16 , 15) 
xt  Z  xf 

(51p) 

^uf  “ 

[Hehf 

(5  lq) 

For 

carbon,  set  X(C)  =  0,  P(C)  =  1. 

tc]uf  = 

[C]xf  +  [C0]xfG(C,16)  +  [C02]  H(15, 16, C)  . 

(51r) 

To  maximize  the  ionization  at  early  time,  as  suggested  by  W. 
Knapp,  we  re-establish  the  molecular  ions  produced  and  dissociatively 
recombined  by  the  Group-X  and  -U  depositions  (and  allow  the  molecular 


ion  decay  to  be  predicted  by  the  later-time  chemistry  module).  First, 
record  the  densities  of  the  molecular  ions  N2,  N0+,  and  C>2  which  under¬ 
went  dissociative  recombination  during  Group-U  deposition. 


[Nj]  =  a,{ [N2]  -  [N2]  } 

u  xf  uf 


(52a) 


[N0+]u  =  [N0]xf{l  -  P (7) } 


(52b) 


[02]  32([02(X)]x£  -  [02  (^)  -^uf^  ~  [®2^^^uf^  "  (52c) 


Next,  add  in  the  corresponding  ions  from  the  Group-X  deposition. 


[M+]  = 

1  d 

[Nt]  = 

1  d 

^[N+j  +  [n2] 

U 

(53a) 

[Mt]  = 
z  d 

[N0+]d  = 

62[no+]  +  [N0+]u 

(53b) 

1 — 1 

1 _ 1 

L 

III 

[ot]  = 

61C°2]  +  62[02]  +  [02]  * 

(53c) 

We  must  subtract  the  products  of  dissociative  recombination 
from  the  appropriate  species  densities,  with  the  constraint  that  N  and 
0  atom  densities  remain  non-negative.  If  limiting  is  required,  main¬ 
tain  the  fractional  contributions  and  scale  appropriately. 

The  total  density  of  these  molecular  ions  is 

[M+]d  =  £  [M it]  (54) 

a  i=l  1  d 

and  the  fractional  contributions  by  each  of  the  ions  is 

ft  =  [Mt]  /[M+]d  ,  i-  1,3  .  (55) 

d 

For  the  dissociative  recombination  reactions  for  Nt,  N0+,  and  02  (see 
Table  4-6),  let  Rp  R2,  and  R-j  be  the  number  of  N(^S),  N(^D),  and  0 
atoms  produced  per  dissociative  recombination  with  fractional 


(56a) 


contributions  by  N^,  N0+,  and  C>2  of 

R]_  =  ^2^2 

R2  =  2(1-B1)f1  +  (l-e2) f 2 

R3  =  f2  +  2f3  . 

For  a  total  density  of  dissociative 
the  densities  of  the  product  atomic 

sid  -  t«Vi  ■ 

The  tentative  final  values  for  thes< 
^iuf  ~  ^iuv  -  ^id  ’ 


f  d ,  f2,  and  f^.  Then 

(56b) 

(56c) 

recombination  equal  to  [M+]d, 
species  are 

1=1,3  .  (57) 

atomic-species  densities  are 

1=1,3  (58) 


where  S^uv  is  given  by  Equation  (51c) ,  (51e) ,  and  (51d)  for 
i  =  1,2,3,  respectively.  If  >  1  for  each  of  the  three  atomic 

species,  no  revision  is  required;  otherwise,  a  new  value  of  [M+] ^ 
is  determined: 


[M+] 


d 


(<Siuv  -  1)/R1 
Min< 

(0.5  S.  /R. 

V  iuv'  i 


S.  >1 
xuv 


S.  <1 

1UV  — 


(59) 


Then,  revised  values  are  computed  for  S^d>  and 

[M*]  =  [M+]df.  ,  1=1,3  .  (60) 

d 

The  final  molecular- ion  densities  are 

[Mt]  =  [Mt]  +  [M+]  ,  1=1,3  .  (61) 

1  uf  1  xf  1  d 

There  are  two  equivalent  ways  of  obtaining  the  atomic- ion 
densities.  One  is  to  solve  the  differential  equations  and  the  other 
is  to  appropriately  difference  the  neutral  species  densities.  The 
results  of  the  two  methods  can  be  reconciled  by  adding  and  subtracting 
appropriate  terms,  followed  by  regrouping.  Here  we  present  the 
results  from  the  second  method. 


[n\£+  [N(4S)]xf  -  [N(4S)]uv 

+  [N(2D)]xf  -  [N(2D)]uv  +  [N(2P)]xf  -  [N(2P)]uf 


+  [NO]  ,  -  [NO]  .p  +  [N02]  -  [NO-] 

xr  u£  xf  Z  uf 


+  2{[N2]  -  [N2]  } 

xf  uf 


[0+]uf  -  [0+]xf  +  [0]xf  -  [0]uv  +  [N0]xf  -  [NO]uf 

+  [OH]xf  -  [0H]uf  +  [C0]xf  -  [CO]uf  +  [H20] 

-  [H20]  +  2{[02(X>]  f  -  [02(X)]uf  +  [02(a>]  f 

uf 

-  [02(a)]uf  +  [N02]xf  -  [N02]uf  +  [H02]xf  -  [H02]u£ 

+  [CO,]  -  [CO,]  )  +  3{[0,]  -  [0,]  )  .  (63) 

xf  uf  xf  uf 

The  final  electron  density  is 

W„£  -  [“+]uf  +  t°+hf  +  [<3uf  ■  <64> 


There  remains  the  task  of  devising  a  truly  satisfactory  pro¬ 
cedure  for  computing  a  pressure  increment  in  addition  to  that  computed 
in  Subroutine  PHEAT.  The  treatment  in  Subroutine  PHEAT  is  incomplete 
because  the  Group-U  subgroups  are  deposited  independently  of  each  other 
across  an  appreciable  path  length  without  accounting  for  either  species 
depletion  by  other  subgroups  or  additional  absorbers  resulting  from  the 
(instantaneous)  dissociative  recombination  that  is  included  in  Sub¬ 
routine  PCHEM.  (Note,  however,  that  the  partial  heating  computed  in 
Subroutine  PHEAT  does  serve  the  purpose  of  defining  effective  fluences 
for  the  Group-U  subgroups  that  are  used  in  Subroutine  PCHEM • )  Any 
additional  pressure  increment  computed  in  Subroutine  PCHEM  must  be  con¬ 
sistent  with  (1)  that  computed  in  Subroutine  PHEAT  and  with  (2)  the 
return  of  pressure  computed  by  the  later- time  chemistry  in  CHEMEF.  The 
problem  is  further  complicated  by  the  re-establishment  of  the  molecular 


72 


ions  produced  by  the  Group-X  and  Group-U  depositions.  To  compensate 
for  (1)  removing  a  number  of  N  and  0  atoms  (whose  sum  is  2[M+]d>  con¬ 
sisting  of  SldN(4S)  atoms,  S2dN(2D)  atoms,  and  S3d  0  atoms)  at  tempera- 


ture  Tnew  and  a  number  [M  ]d  of  electrons  at 
returning  energies  (in  eV)  of 

temperature  Tgf,  and  (3) 

e(n2) 

=  15.58  -  9.76  -  2(1-8^2.37  = 

3.45 

(65a) 

e(02) 

=12.06-5.12 

6.94 

(65b) 

E  (NO"*") 

9.27  -  6.51  -  (l-e2)  2.37  = 

0.9825 

(65c) 

+  +  + 

each  time  a  re-established  C^,  and  NO  ion  dissociatively  recom¬ 

bines,  we  should  subtract  an  amount  of  thermal  energy  equal  to 

[M+]d 

{ 1 . 5k (Tnew  +  Tef  -  2Told)  +  Ed} 

eV/cm3  , 

(66) 

where 

-  +  E  (O^)  f  2  +  E(N0+)  f . 

• 

(67) 

However,  as 
given  by 

an  interim  approximation  we  shall 

account  for  only 

AP 

AP 

=  -[M+]d  Ed  1.6  >'  10" 12 

3 

erg/cm 

(68) 

pi 

=  P  .  .  +  AP  . 
old 

(69) 

To  ensure  that  the  heavy-particle  pressure  will  always  be 
positive,  we  have  added  the  electron  pressure  (computed  below)  to  the 
previously-found  value  for  the  total  pressure.  This  interim  measure 

needs  to  be  reviewed  and  improved. 

To  compute  the  temperature  TEF  taken  to  be  common  to  the 

1  3 

thermal  electrons,  the  N2  vibrational  states,  and  the  0(  D)-to-0(  P) 
population  ratio,  we  use  the  GET-prepared  Subroutine  TEXK  [Vol.  11] 
which  determines  srch  a  temperature  when  given  [ N2 3  >  [0] ,  [e],  and  the 
energy  per  cubic  centimeter  (ET)  to  be  divided  among  the  three  modes. 
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As  an  interim  measure  we  compute  ET  from 


ET  -  2  *  5M*,U  +  ^HP  *  Epe 


where  6[e]x  is  the  incremental  electron  density  produced  by  the 
Group-X  and  -U  depositions,  the  factor  2  is  an  interim  assumption  for 
the  energy  in  the  three  modes  per  newly-formed  electron  in  the  Group-X 
and  -U  depositions,  is  the  incremental  electron  density  produced 

by  the  heavy-particle  depositions,  and  Epe  is  the  sum  of  the  electron 
thermal  and  0(^D)  excitation  energies  per  newly-formed  electron  in  the 
heavy-particle  deposition,  obtained  from  Subroutine  HPCHEM. 

Sum  the  final  densities  of  N,  0,  C,  and  H  nuclei,  followed 
by  finding  the  fractional  increases  in  the  final  values  relative  to  the 
initial  values,  as  measures  of  conservation  errors  in  treating  the 
Group-U  deposition. 

The  foregoing  describes  the  combined  treatment  of  target 
cells  that  are  Case  1  or  Case  3.  For  Case  1,  the  treatment  is  complete, 
but  for  Case  3,  the  final  cell  properties  are  determined  by  volume¬ 
weighting  the  two  portions,  one  of  which  is  inside  the  B-edge  and  the 
other  is  outside  the  B-edge. 

For  Case  2,  the  entire  cell  is  inside  the  B-edge  and  the 
only  remaining  species  are  assumed  to  be  N+,  0+,  and  e.  The  tempera¬ 
ture  of  both  electrons  and  heavy  particles  is  taken  to  be  1  eV. 

The  final  step  in  Subroutine  PCHEM  is  to  store  the  cell 
quantities  in  temporary  storage  BUF1  until  the  loop  over  the  entire 
column  of  target  cells  is  completed,  after  which  the  entire  column  of 
cell  quantities  is  transferred  to  scratch  storage  in  large-core  memory. 
After  the  energy  deposition  is  completed  for  the  entire  grid,  cell 
quantities  in  the  scratch  area  are  transferred  to  Time-Slot  2. 
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APPENDIX  A 


SPECIES  FROM  X-RAY  DEPOSITION 


A-l  INTRODUCTION 

A  procedure  for  calculating  the  species  formed  following  x- 
ray  deposition  is  presented.  This  procedure  is  generally  an  improve¬ 
ment  over  that  in  ROSCOE-Radar  of  assuming  ions  (of  neutral  species) 
are  formed  in  proportion  to  the  densities  of  neutral  species  [HS-75b, 
p.  108].  Here  we  take  account  of,  not  only  collisional  ionization, 
but  also  dissociative  ionization,  collisional  dissociation,  colli¬ 
sional  excitation,  and  the  dependence  of  ion  production  on  arbitrary 
mixes  of  the  neutral  species.  In  addition,  a  more  accurate  extimate 
of  the  CO2  loss  resulting  from  x-ray  deposition  is  given. 

A- 2  PRINCIPAL  SPECIES  FORMATION 


In  this  section  we  will  initially  consider  ion  production 
from  N2,  O2,  and  0  as  absorbers;  later,  we  include  N  as  an  absorber 
and,  in  some  parts  of  the  calculation,  NO.  (The  full  treatment  of  NO 
and  eight  other  minor  species  is  discussed  in  Section  4.)  Later  in 
this  section  we  consider  the  neutral  species  from  non-ionizing  colli¬ 
sions  of  the  slowing-down  electrons. 

To  predict  the  ion  species  resulting  from  x-ray  deposition, 
we  modify  the  relations  between  ion  species  and  fractional  neutral 
density  previously  developed  [MS-76a]  for  photon  absorption.  These 
relations  are  of  the  form  [cf.  MS-76a,  Equation  (16)] 


a . 


6, 


<s  ;i>  ■  aifil  +  Vi1  VEP 


(Al) 


where  (S+; i)  is  the  number  of  ions  formed  from  absorber  i  per  photon 
absorbed  in  the  neutral-particle  mix  of  N£,  C^,  and  0.  For  molecular 
absorbers,  (S+;i)  is  the  sum  of  the  resulting  molecular  and  atomic 
ions,  assumed  to  be  singly  charged.  The  quantities  a,  b,  a,  3,  and 
are  constants  [MS-76a,  Table  6,  p.  36],  Ep  is  the  energy  of  the 
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absorbed  photon,  and  is  the  fractional  density  of  absorber  i  among 
the  neutral  species.  In  MS-76a,  the  values  of  the  constants  were 
developed  for  the  neutral  species  N2.  O2.  and  0  on  the  basis  of  photon 
deposition  results  computed  for  ambient  atmospheres  with  compositions 
given  by  the  CIRA-1972  mean  atmosphere  [Cl- 72]  corresponding  to  the 
range  of  altitudes  from  60  to  250  km  and  for  photon  energies  of 
33.67  eV  and  larger.  These  constants  are  reproduced  here  in  Table  A-l. 

In  using  Equation  (Al)  for  present  purposes,  we  need  to 
address  three  problems:  (a)  extending  the  use  of  Equation  (Al)  for 
photon  energies  greater  than  281  eV,  (b)  extending  the  use  of  Equation 
(Al)  for  chemically-altered  atmospheres,  and  (c)  incorporating  N  atoms 
as  absorbers. 

A- 2 . 1  Ion  Formation  Following  High-Energy  Photon  Absorption  in  an 

Ambient  Atmosphere 

The  derivation  of  Equation  (Al)  and  the  constants  of  Table 
A-l  were  based  on  calculations  in  which  the  maximum  photon  energy  was 
281  eV  [MS-76a].  To  calculate  and  verify  the  accuracy  of  the  ions  per 
absorbed  photon  at  higher  photon  energies,  we  take  two  steps.  Firstly, 
we  compute  the  ions  per  absorbed  photon  using  Equation  (Al)  but  exclud¬ 
ing  the  ions  formed  in  the  x-ray  absorption  event.  Thus,  the  counted 
ions  per  absorbed  photon  result  from  the  slowing-down  of  a  photoelec¬ 
tron  formed  in  the  x-ray  absorption  event.  Such  ions  per  absorbed 
photon  are  shown  to  be  accurate  by  comparison  with  detailed  calculations 
[MS- 75,  MS- 76a]  on  the  slowing-down  of  high-energy  electrons.  Secondly, 
we  modify  Equation  (Al)  to  include  the  results  of  the  x-ray  absorption 
event  for  photon  energies  above  281  eV. 

Equation  (Al)  was  derived  from  calculations  of  the  absorption 
of  photons  with  energies  below  the  K-absorption  edge  of  both  N  and  0 
atoms.  According  to  the  x-ray  deposition  model  previously  formulated 
[MS-75,  Appendix  B] ,  the  absorption  event  in  this  case  resulted  in  one 
electron  and  one  ion.  Therefore,  to  use  Equation  (Al)  to  describe  only 
the  effect  of  the  photoelectron  slowing-down,  it  must  be  modified  as 
follows : 


=  (S+;i) 


fi 


A-  2 


1 


1 


d 


(S+;iK 


(A2) 


Table  A-l.  Constants  used  in  Equation  (Al)  relating 
absorber  and  species  formed  per  absorbed 
photon  [MS-76a,  Table  6]. 


Absorber  i- 

a . 

X 

bi 

ai 

E;.  eV 

N2 

1.06 

0.0266 

0.872 

0.649 

33.67 

°2 

1.26 

0.0238 

0.927 

0.821 

33.67 

0 

1.20 

0.0350 

1.09 

1.30 

33.67 

—  Later,  in  Section  A-2.2.3  when  we  add  nitrogen  atoms  as 
absorbers,  we  shall  postulate  that  the  constants  for  N 
are  the  same  as  those  for  0. 


Here,  (S+;i)^  is  the  number  of  ions  formed  from  absorber  i  per  photo- 
electron  and  (S+;i)  and  are  as  defined  above.  The  use  of  f^  as 
opposed  to  the  accurate  quantity,  0-|_f  ^  (where  0 is  the  photo¬ 

absorption  cross-section) ,  is  an  approximation  which  introduces  an 
error  of  less  than  5  percent  in  the  case  of  the  slowing  down  of  a  900- 
eV  photoelectron.  Equation  (A2)  will  now  be  tested  for  higher-energy 
(>281  eV)  photons  against  the  detailed  calculations  [MS-75,  MS-76a]. 

To  use  Equation  (A2)  in  the  comparison  with  the  detailed 
calculations,  the  contributions  of  atomic  and  molecular  ions  to 
Equation  (A2)  must  be  separated.  This  separation  is  accomplished  with 
the  following  relationships: 


and 


Rn+  =  (N+)/[(N+)  +  (N+)] 


=  0.184  +  5.8  x  10"5  (z-60) 

%  0.192  for  z  =  200  km 


(A3) 


R0+  =  (0+)/[(0+)  +  (0+)]  =  0.302 


(A4) 


A-  3 


where  z  is  the  altitude  (km) .  Both  of  these  equations  are  based  on 
calculations  of  electron  slowing-down  [MS-76a,  Tables  10a-10d].  Equa¬ 
tion  (A3)  gives  (within  about  2  percent)  the  ratio  of  the  N  ions  to 
the  total  ions  derived  from  ^  for  incident  electron  energies  in  the 
range  from  230  to  900  eV  and  for  atmospheric  compositions  corresponding 
to  the  range  of  altitudes  from  60  to  250  km.  (See  Tables  10a,  10b,  10c, 
and  lOd  of  MS-76a.)  Equation  (A4)  gives  the  ratio  of  the  cf*~  ions  to  the 
total  ions  derived  from  C^,  averaged  over  the  range  of  incident  electron 
energies  from  230  to  900  eV. 

In  making  the  comparison  between  the  predictions  of  Equation 
(A2)  and  the  detailed  calculations  [MS-76a],  we  consider  a  range  of 
energies  of  the  x-ray  photon  and  the  incident  electron.  For  the  de¬ 
tailed  calculations  [MS-76a],  the  selected  range  for  the  incident  elec¬ 
trons  is  from  230  to  900  eV ;  for  the  calculations  based  on  Equation 
(A2) ,  the  selected  range  for  the  incident  photons  is  from  650  to  1400 
eV.  For  this  latter  range,  the  photon  energy  will  be  above  the  K- 
absorption-edge  energy  for  N  and  0  atoms  and  will,  generally,  yield 
photoelectrons  having  energies  within  the  range  used  in  the  detailed 
calculations.  However,  the  ranges  of  energies  are  not  explicitly  con¬ 
sidered  in  the  comparison  since  the  quantity  compared,  the  ion  species 
per  ion  pair,  does  not  vary  significantly  over  the  energy  range  speci¬ 
fied  for  either  the  detailed  calculations  or  those  based  on  Equation 
(A2) .  The  variation  is  generally  not  greater  than  5  percent  in  the 
former  calculations  and  generally  not  greater  than  1  percent  in  the 
latter  calculations . 

The  results  of  the  calculations  are  given  in  Table  A-2A  in 
terms  of  the  ion  species  per  ion  pair.  This  quantity  is  calculated 
from  Equation  (A2)  by  forming  the  ratio 

(S+;i)'p  =  <S+;i)7  £  (S+;J)'  .  (A5) 

i 

The  entries  are  averages  over  the  energy  ranges  cited  above  for  the 
detailed  calculations  and  those  based  on  Equation  (A2) .  For  all  com¬ 
parisons,  the  relative  error  is  not  greater  than  about  16  percent;  for 
the  majority  of  cases,  the  relative  error  is  less  than  9  percent. 
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Table  A-2.  Comparison  of  predictions  using  Equation  (A2)  with 
earlier,  detailed  calculations. 


A. 

Comparison 

of  ion  species  per 

60  km 

ion  pair  (columns  2  and 

110  km 

3) 

Species 

Eq.  (A2) 

MS- 76a- 

7o  Error 

Eq.  (A2) 

MS-76a- 

o 

/ 

i  Error 

-j 

0.632 

0.629 

0.5 

0.611 

0.624 

-2.1 

N+ 

0.142 

0.142 

0.0 

0.141 

0.143 

-1.4 

°2 

0.158 

0.156 

1.3 

0.102 

0.091 

12.1 

0+ 

0.068 

0.073 

-6.8 

0.147 

0.143 

2.8 

145  km 

250  km 

Species 

Eq.  (A2) 

MS- 76a- 

70  Error 

Eq.  (A2) 

MS- 76a- 

2 

o  Error 

N2 

0.459 

0.502 

-8.6 

0.198 

0.217 

-8.8 

N+ 

0.107 

0.117 

-8.5 

0.048 

0.052 

-7.7 

°i 

0.053 

0.051 

3.9 

0.012 

0.013 

-7.7 

(T 

0.382 

0.330 

15.8 

0.741 

0.719 

3.1 

B.  Comparison  of  electron 
electrons) 

volts  per  ion  pair 

(for  900-eV 

Altitude 

(km) 

60 

110 

14  5 

250 

Eip (pe) -  34.8 

35.5 

32.8 

28.1 

34.1 

32 . 5 

29.2 

—  Results  for  60- ,  110- ,  145- ,  and  250-km  altitude  are  obtained  from 
an  average  over  the  range  230  £  U(eV)  <_  900  of  the  results  from 
Tables  10a,  10b,  10c, and  lOd,  respectively,  in  MS-76a. 

-7.  error  =  {[Eq.  (A2))  -  (MS- 76a)  ]/ (MS- 76a) }  *  100- 
c  + 

—  Calculated  by  using  E.  (pe)  =  900/(E(S  ;i)-l);  also  see  Equation 

(19)  of  MS-76a.  ip  i 

—  From  Table  11a  of  MS- 76a  for  U  =  900  eV. 


In  Table  A-2B,  values  of  the  electron  volts  per  ion  pair, 

E.  ,  are  given  for  the  detailed  calculations  (E.  (e))  and  those  based  on 

ip  JLp 

Equation  (A2)  (E.^  (pe)  )  .  For  this  comparison,  both  the  photoelectron 
in  the  x-ray  absorption  case  and  the  incident  electron  in  the  detailed 
calculation  had  an  energy  of  900  eV.  The  agreement  is  adequate.  There 
is  some  uneveness  in  the  E^p(pe)  values  as  a  function  of  altitude  and 
this  is  attributed  to  the  compromises  in  the  fitting  of  the  parameters 

» 

in  Equation  (Al)  to  the  data  from  the  detailed  calculations  [MS-76a] 
on  photon  absorption. 

*  From  the  above  calculations  and  comparisons,  we  find  Equa¬ 
tion  (A2)  to  be  satisfactory  in  computing  the  ion  densities  resulting 
from  photoelectron  slowing-down  in  ambient  atmospheres  following 

% 

absorption  of  x-rays  with  energies  of  650  eV  and  larger. 

* 

4.  The  ion  species  per  ion  pair  in  the  given  neutral-particle 

mix,  (S+).  ,  resulting  from  absorption  of  x-rays  with  energies  above 


i  650  eV,  are  given  by  the  following  equations: 

r  (nJ)±  =  [(S+;N2)  -  2fw  JRn+/IP  (A6 

i  P  2  2 

(N+)ip  =  [(S+;N2)RN+  +  2fN2RN+]/IP  (A7 ) 

^  <0+>  =  [(S+;02)  -  2fQ2]R0+/IP  (A8) 

(0+).p  =  [(S+;02)R0+  +  2fQ2R0+  +  (S+;0)]/IP  (A9) 

where 

IP  =  L  (S+;i)  (A10) 

i 

rn+  =  1  -  V  (Alla) 

Rq£  1  =  R0+  .  (Allb) 
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In  Equations  (A6)  through  (A9) ,  the  factor  2  arises  from  the  assumption 
that  in  the  photoabsorption  event  two  atomic  ions  and  two  photoelectrons 
are  formed  from  each  molecular  absorber  per  absorbed  photon. 


In  Table  A-3  we  compare  the  results  obtained  from  Equations 
(A6)  through  (A9)  with  (a)  the  detailed  calculations  modified  to  account 
for  an  x-ray  absorption  event  (by  adding  the  ion  corresponding  to  the 
photoelectron)  and  with  (b)  the  ion  species  obtained  from  the  prescrip¬ 
tion  used  in  ROSCOE-Radar  [HS-75b].  The  new  method  of  calculating  the 
ion  species  per  ion  pair,  Equations  (A6)  through  (A9) ,  is  shown  in  Table 
A-3  to  be  a  generally  better  approximation  to  the  detailed  calculations 
than  the  older  presciption  [HS-75b].  Table  A-3  applies  to  absorptions 
in  ambient  atmospheres  (in  contrast  to  a  chemically-disturbed  atmosphere 
considered  in  Section  A-2.2). 


Equations  (A6)  through  (A9)  depend  on  the  photon  energy,  E  , 
through  Equation  (Al) .  However,  for  the  number  of  ions  per  ion  pair 
the  dependence  on  photon  energy  is  sufficiently  weak,  provided  E^  is 
greater  than  650  eV,  to  justify  eliminating  this  dependence.  If  E^  in 
Equation  (Al)  is  fixed  at  1400  eV,  then,  for  photon  energies  greater 
than  650  eV,  the  error  in  calculating  the  ion  species  from  the 
specialized  equation  will  not  exceed  about  10  percent  compared  with 
accounting  for  the  energy  dependence.  Making  this  simplification  leads 
to  the  following  equation: 


(S+;i).  =  (a.f.^  +  b'f^i)/IP 

’  ip  'll  1  i  ' ' 


IP 


(A12) 


where 


b:  =  1366  b± 


(A13) 


A-2.2  Ion  Formation  in  Chemically-Disturbed  Atmospheres 

To  calculate  the  ion  species  per  ion  pair  in  chemically- 
disturbed  atmospheres,  consideration  must  be  given  to  the  altered  com¬ 
position  as  compared  with  ambient  atmospheres. 


Table  A-3.  Comparison  of  ion  species  per  ion  pair  calculated 
according  to  (a)  the  new  method,  (b)  the  detailed 
method  [MS-76a],  and  the  older  ROSCOE-Radar  method 
[HS-75b].£ 


z,km 

Method 

(N't)  . 

v  2  '  ip 

7o  E- 

(N+)  . 

V  7  ip 

%  E 

(Ot)  . 

v  2'  ip 

%  E 

<0+>ip 

7o  E 

60 

new 

0 

.618 

1. 

.3 

0, 

156 

-2, 

,5 

0 

.154 

0, 

.  7 

0, 

.072 

-7, 

.7 

detailed 

0 

.610 

0. 

160 

0, 

.153 

0. 

.078 

old 

0 

.790 

29. 

.5 

0. 

000 

-100 

0. 

.210 

37, 

.3 

3 

<10-6 

-100 

110 

new 

0 

.597 

-2, 

.0 

0. 

155 

-3. 

.  7 

0 

.099 

10 

.0 

0. 

.148 

5, 

.0 

detailed 

0 

.609 

0. 

161 

0. 

.090 

0. 

.141 

old 

0 

.730 

19. 

.9 

0. 

000 

-100 

0 

.120 

33. 

.3 

0. 

.150 

6. 

.4 

145 

new 

0 

.450 

-9. 

5 

0. 

117 

-10. 

7 

0. 

.052 

^0 . 

,0 

0. 

,381 

19. 

.4 

detailed 

0 

.497 

0. 

131 

0. 

.052 

0. 

.319 

old 

0 

.526 

7. 

.3 

0. 

000 

-100 

0. 

.059 

15. 

.7 

0. 

,415 

30. 

.1 

250 

new 

0. 

.194 

-12. 

,2 

0. 

051 

-12. 

0 

0. 

,012 

-8. 

,3 

0. 

,743 

4. 

,8 

detailed 

0. 

.221 

0. 

058 

0. 

013 

0. 

709 

old 

0. 

.183 

-17. 

2 

0. 

000 

-100 

0. 

012 

-8. 

3 

0. 

805 

13. 

5 

a 

-  The  ion  species  per  ion  pair  derived  from  the  new  method  apply  strictly 
to  a  1400-eV  x-ray  photon;  the  results  presented  here  for  the  detailed 
method  have  been  averaged  over  those  for  450-eV  and  900-eV  photo¬ 
electrons  and  take  account  of  the  x-ray  absorption  event  (by  adding  the 
ion  corresponding  to  the  photoelectron) .  These  two  calculations  are 
regarded  as  comparable  since,  for  the  photoelectron  slowing-down ,  the 
resulting  ion  species  per  photoelectron  are  insensitive  to  the  photo¬ 
electron  energy,  at  least  for  energies  of  450  eV  and  larger. 

-  70  E  =  percent  error  =  {[(new  or  old)  -  (detailed)  ]/ (detailed}  *  100. 
Note  that  with  this  definition  a  negative  error  is  formally  limited 
to  -1007o. 


A-2.2.1  Effect  of  Ambient  Electrons 


The  effect  of  ion  formation,  due  to  thermal  electron  densities 
significantly  greater  than  ambient  values,  has  been  determined  [MS-76a, 
pp.  61,63,64]  for  high-energy  electron  deposition.  For  high-energy 
electrons  (450  and  900  eV) ,  the  ion  species  per  ion  pair  were  changed 
by  less  than  about  10  percent  for  electron  densities  10^  times  the 
ambient  value  at  altitudes  in  the  range  from  100  to  250  km  and  about 
1011  times  the  ambient  value  at  60  km.  At  the  upper  limit  to  the  elec¬ 
tron  density  in  these  studies,  the  steady-state  approximation  for  the 
electron  energy  distribution,  used  in  the  calculations,  became  invalid 
because  of  the  loss  of  non-thermal  electrons  via  dissociative  recom¬ 
bination.  Maintaining  the  validity  of  the  steady-state  approximation 
limited  the  ratio  of  thermal  electrons  to  neutrals  at  60,  110,  and  145 
km  to  1,5  x  10~\  5  x  10"^,  and  3  x  10"^,  respectively;  at  250  km,  an 
electron  density  of  10^  times  the  ambient  value  corresponded  to  an 
average  of  2.6  free  electrons  per  neutral  particle.  However,  at  250 
km,  where  the  ratio  of  thermal  electrons  to  neutrals  was  not  signifi¬ 
cantly  limited,  the  change  in  ion  species  per  ion  pair  was  of  the  order 
of  10  percent  or  less;  thus,  it  is  reasonable  to  assume  that  species 
formation  by  electron  or  photon  deposition  is  not  strongly  affected  by 
an  increased  electron  density  even  up  to  72  percent  ionization  at  all 
altitudes  considered . 


A-2.2.2  Effect  of  Ambient  Ions 


The  effect  on  ion  formation,  due  to  ion  densities  signifi¬ 
cantly  greater  than  ambient  values,  has  not  been  determined.  A  simple 
prescription  to  partially  account  for  their  effect  is  developed  here. 

Consider  the  disturbed  atmosphere  to  consist  of  an  ion  gas 
(of  singly-ionized  particles)  and  a  neutral  gas.  The  fraction  of  the 
x-ray  energy  deposited  in  the  neutral  gas  will  be  modeled  by  assuming 
that  the  probability  that  an  electron  expends  energy  on  a  target  parti¬ 
cle  is  proportional  to  the  number  density  of  the  particle,  inversely 
proportional  to  its  ionization  energy  I  (or  dissociation  energy,  for 
molecular  ions) ,  and  inversely  proportional  to  its  effective  mean 
excitation  energy  E  defined  by  the  expression 


g 


1 

E/I 
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(A14) 


where  E^  is  the  energy  expended  (eV)  in  forming  an  ion  pair  within  the 
given  particle  gas.  For  neutral  and  ion  gases  consisting  of  atomic  and 
molecular  species,  the  portion  of  the  total  energy  deposited  per  unit 
volume,  E^,  which  is  assigned  to  the  neutral  gas,  E^n,  is  given  by 

V  T 

—  n 

(A15) 


Edn'Ed 


n 


n 


ETn  +?T 

n  1 


where 


_  [n] 


n 


n 


a+gn> 


n  =  N, 


NO,  N,  0 


(A16) 


T.  =  (1+gp  ,  i  =  N2,  0+,  N0+,  N+,  0+  .  (A17) 

Values  of  the  parameters  entering  Equations  (A14) ,  (A15) , 

(A16) ,  and  (A17)  are  given  in  Table  A-4.  The  value  of  g  for  the  atomic 

species  is  taken  to  be  unity;  for  0,  this  value  is  based  on  the  fact 

that  the  energy  per  ion  pair  in  an  0  gas  is  about  twice  the  ionization 

energy  [MS-76a]  and  for  N,  0+,  and  N+,  on  the  assumption  of  similar 

values.  The  value  of  g  for  neutral  molecules  is  taken  to  be  4/3  on 

the  basis  that  E.  /I  is  2.24  for  N0  and  2.48  for  0o ;  for  molecular  ions, 
ip  +  +2  +2 

g  is  taken  to  be  3/4  for  N2  and  O2  and  4/3  for  NO  under  the  assumption 
that  E  is  about  the  same  as  for  neutral  molecules.  An  improved  basis 
for  assigning  these  values  would  certainly  be  desirable. 

A-2.2.3  Addition  of  Nitrogen  Atoms  as  Absorbers 

In  chemically-disturbed  atmospheres,  the  density  of  nitro¬ 
gen  atoms  may  increase  to  represent  a  significant  set  of  absorbers  of 
the  x-ray  photons.  To  account  for  this  situation,  we  add  nitrogen 
atoms  to  the  set  of  absorbers  by  modifying  Equations  (A6)  through  (A9) . 

Since  calculations  of  the  absorption  and  slowing-down  process 
have  not  been  performed  for  atmospheres  containing  nitrogen  atoms,  we 
assume  that  nitrogen  atoms  may  be  treated  analogously  to  oxygen  atoms 
for  the  purpose  of  calculating  the  ion  species  per  ion  pair  from  Equa¬ 
tion  (A12) .  In  this  case  Equations  (A6) ,  (A8) ,  and  (A9)  remain 


Table  A-4.  Parameters  used  in  the  expressions  for 

partitioning  the  deposited  energy  between 
the  neutral  and  ionized  gases. 


Species 

I-,  eV 

Eip-.  eV/ip 

n2 

15.6 

4/3 

35 

°2 

12.1 

4/3 

30 

NO 

9.27 

4/3 

25 

N 

14.5 

1 

27 

0 

13.6 

1 

27 

»2 

8.7 

3/4 

°2 

6.7 

3/4 

N0+ 

10.8 

4/3 

N+ 

29.6 

1 

0+ 

35.1 

1 

2i 

—  I  is  the  ionization  energy  except  for  the  molecular 
ions  for  which  it  is  the  dissociation  energy.  Values 
taken  from  BB-74a,  pp.  84,85. 

—  Values  of  the  energy  per  ion  pair  are  nominal  selections 
from  those  computed  in  MS- 76a  for  No,  Oo,  and  0.  The 
value  for  N  is  assigned  to  be  that  for  0.  The  value  for 
NO  is  a  postulated  value. 

p 

—  See  text  which  discusses  the  assignment  of  these  values. 


unchanged,  but  Equation  (A7)  becomes  analogous  to  Equation  (A9) ,  i.e., 


(N+)ip  =  [(S+;N2)Rn+  +  2fN2RN+  +  (S+;N)]/IP  (A7a) 

where  the  constants  used  in  Equation  (Al)  to  evaluate  (S+;0)  will  be 
taken  to  be  those  used  to  evaluate  (S+jN),  and  the  sum  in  Equation 
(A10)  is  now  over  the  species  i  =  l$2>  ,  0,  and  N. 

A-2.2.4  Mean  Energy  Per  Ion  Pair 

The  energy  per  ion  pair  used  in  ROSCOE-Radar  [HS-75b,  p.  108] 
to  compute  the  number  of  ion  pairs  from  the  energy  deposited  is  35  eV/ip. 
At  high  altitudes  where  atomic  oxygen  dominates  the  atmospheric  com¬ 
position,  the  use  of  this  value  can  introduce  an  error  of  about  20  per¬ 
cent.  A  more  accurate  value  for  the  energy  per  ion  pair  in  the  neutral 
gas  may  be  obtained  from  the  prescription 


fsEip<s> 


(A18) 


where  values  of  E^(s)  are  given  in  Table  A-4  for  s  =  N2,  C^,  NO,  N, 
and  0.  In  atmospheres  without  NO  and  N,  Equation  (A18)  is  accurate  to 
within  about  3  percent. 


A- 2. 3  Neutral  Species  Formation  from  Non-Ionizing  Collisions 

MS-76a  (p.  36)  notes  that  relations  similar  to  Equation  (Al) 
could  be  developed  for  neutral  species,  but  that  no  attempt  was  made  to 
do  so.  In  the  absence  of  such  relations,  we  shall  resort  to  simpler 
procedures.  We  are  mainly  interested  in  the  atoms  formed  when  slowing- 
down  electrons  collisionally  dissociate  ^  and  C^.  We  shall  also  con¬ 
sider  the  excitation  of  02(a)  since  this  species  is  carried  in  the  late¬ 
time  chemistry  used  in  the  high-altitude  grid. 

Let  Fj  denote  the  ratio  of  collisional  dissociation  to  colli- 
sional  ionization,  i.e., 


Fd(N2) 

N2  collisional  dissociation 

N2  collisional  ionization 

(A19a) 

Fd(°2) 

C>2  collisional  dissociation 

02  collisional  ionization 

(A19b) 
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From  the  data  in  MS-76a  (Tables  10a  through  lOd)  one  can  infer  that, 
for  900-eV  electrons, 


Fd(N2)  v  0.22  to  0.31, 

a  range  which  depends  on  the  way  predissociation  states  are  treated. 
This  inferred  value  of  Fd(N2)  was  obtained  by  using  the  distribution 
of  species  per  ion  pair  to  evaluate  the  expression 


0.5[£N(s)  -  N+] 


W  - 


4  2  2 

where  s  =  S,  ZD,  ZP. 


(A20) 


F.  Gilmore  has  begun  a  review  of  the  question  of  predissoc¬ 
iation  in  N2  and  advocates  more  excited- state  nitrogen  atoms  per  ion 
pair  than  indicated  in  MS- 76a.  Gilmore  [DNA  IR  Data  Review  Meeting  at 
AFGL,  11  April  1979]  suggested  0.9  N*  atoms  per  ion  pair.  In  a 
private  communication  on  23  April  1979,  Gilmore  stated  his  review  was 
continuing  but  that  his  best  guess  as  of  the  moment,  for  the  number  of 
nitrogen  atoms  per  ion  pair,  including  dissociative  ionization  as  well 
as  collisional  dissociation,  was 


0.8  NCS)  , 


0.6  N(  D), 


and  0.3  N(ZP)  . 


If  we  use  a  value  of  0.14  IT  per  ion  pair  (obtained  in  MS-75,  MS-76a, 
and  SK-75,  even  though  Gi-74  quotes  0.18)  and  0.63  N2  ions  per  ion 
pair  (as  obtained  in  MS-75,  MS-76a,  and  SK-75;  Gi-74  quotes  0.62), 
then 

Fd(N2)  =  °-5[0._8_.+.  0-6^+  0.3  -  0.14]  =  x  24  ,  (A21) 

larger  than  that  indicated  in  MS- 76a  by  a  factor  of  5.6  to  4.0. 

It  should  be  noted  that  Gilmore's  guessed-values  are  indeed 
preliminary,  and  that  it  seemed  he  was  assuming  0.2  N(^S)  atoms  from 
dissociative  ionization  (in  contrast  to  our  use  of  0.14  N(S>)).  How¬ 
ever,  this  difference  doesn't  change  the  value  of  Fd(N2)  much,  as  seen 
from 
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A..J 


Fd(N2) 


0.5r0.8  -  (0.2-0.14)  +  0.6  +  0.3  -  0.14] 

0.63 


=  1.19. 


(A22) 


Though  it  doesn't  change  the  estimate  of  Fd(N2) ,  we  note  that 
a  possible  assignment  of  dissociating  states  per  ion  pair  is 

0.2N(4S)  +  0.2N+ 

0.3N(4S)  +  0 . 3N(2D) 

0 . 3N(4S)  +  0. 3N(2P) 

0 . 15N(2D)  +  0 . 15N(2D)  . 

We  don't  know  what  Gilmore  had  in  mind. 

We  need  a  prescription  for  assigning  the  relative  abundances 
of  the  nitrogen  atom  states.  Values  inferred  by  averaging  values  from 
MS-76a  [Table  10a,  p.  57]  are 

N(4S)  :  N(2D)  :  N(2P)  =  0.217  :  0.266  :  0.0013 

=  0.448  :  0.549  :  0.003  .  (A23) 

From  Gilmore's  statements  we  have 

N(4S)  :  N(2D)  :  N(2P)  =  0.8  :  0.6  :  0.3 

=  0.471  :  0.353  :  0.176  .  (A24) 

Turning  to  02>  one  can  infer  from  the  data  in  MS-76a  (Tables 
10a  through  lOd)  that  for  900-eV  electrons, 


Fd(02)  =  0.88  .  (A25) 

Based  on  the  900-eV  electron  results  in  Tables  10a,  b,  c,  d 
in  MS-76a,  one  can  show  that  the  following  fit  function  is  satis¬ 
factory: 
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(A26a) 


y  =  {02(a)/0+}ion  pair 
x  =  [0]/[02]  (A26b) 


=|4.8  ,  x  <  0.6528 

^3 . 486  x"3/4  ,  x  >  0.6528  .  (A27) 

A- 3  C02  LOSS  FOLLOWING  X-RAY  DEPOSITION 

The  procedure  in  ROSCOE-Radar  for  treating  the  loss  of  CO2 
following  x-ray  deposition  is  to  assume  the  loss  is  proportional  to 
that  of  N2 .  Here,  a  more  accurate  estimate  is  made  by  determining 
the  relative  loss  of  CC^  on  the  basis  of  electron  flux  spectra  and 
known  electron  energy- loss  functions. 

The  rate  of  loss  of  species  S,  in  the  framework  of  an  energy 
bin  structure  over  the  energy  range,  is  given  by 

-  =  £<Kj)L(j,S)[S]  (A28) 

j 

where  j  is  the  group  number  for  the  energy  bins  [MS- 75,  p.  11],  4>(j)  is 

-  2  - 1 

the  electron  flux  per  unit  energy  in  energy  bin  j  (electrons  cm  sec 

-1  2 
eV  ),  and  L(j,S)  is  the  energy- loss  function  (eV  cm  per  electron)  of 

species  S  in  energy  bin  j.  The  usual  energy- loss  function  has  been 

modified  by  removing  the  contribution  of  processes  not  leading  to  C02 

or  N2  loss.  The  energy-loss  function  values  are  given  in  Table  A-5; 

the  modified  values  are  in  parentheses. 

Values  of  the  species-loss  frequency,  £j<Kj)L(j ,S) ,  are 
given  in  Table  A-6  for  various  altitudes  and  electron  densities,  for 
both  C02  and  N2<  The  electron  flux  values,  <f>(j),  used  in  computing  the 
species-loss  frequencies  were  taken  from  the  calculations  in  MS-76a 
(data  extracted  from  the  microfiche  records  of  these  calculations) . 

The  ratio  of  the  species- loss  frequency  for  C02  to  that  for  N2  can  be 
adequately  represented  by  the  following  expression: 


2>(j)L(j,C02) 

^ -  =  1.92  +  0.7  (1  -  e"1*36*'® 

£  <Kj)L(j,N2) 

j 


F(a) 

(A29) 
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Table  A-5.  Energy- loss  function  values  for  C09  and  N 
up  to  1  keV.  z 


Electron  Energy 
eV 

Loss  Function, 

10"16  cm2  eV 

C02- 

N2b 

11 

%o 

10  (%0)- 

14 

1  (0.9)- 

13  (0.4) 

18 

4 

14  (3.4) 

23 

20 

22  (12) 

29 

34 

31  (21) 

36 

50 

42  (32) 

45 

75 

52  (40) 

57 

85 

63  (50) 

72 

100 

69  (56) 

90 

109 

71  (58) 

110 

120 

71  (59) 

140 

121 

69  (59) 

180 

120 

63  (54) 

230 

118 

57  (49) 

290 

109 

50  (43) 

360 

103 

45  (39) 

450 

99 

39  (34) 

570 

87 

33  (29) 

720 

80 

27  (24) 

900 

70 

23  (20) 

-  SS-72d. 

-  MS-76a,  p.  62 . 


Q 

—  Energy- loss  function  excluding  processes  not 
leading  to  loss  by  target  molecule,  i.e., 
energy- loss  function  for  destruction  of 
target  molecule. 
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(10)  means 


where  the  fractional  ionization,  a,  is  given  by 


e 

t 


a 


t*T 


[e] 

x 


(A30) 


ttl 

[e]  is  the  electron  number  density,  and  [i]  is  the  density  of  the  1 
neutral  species.  Thus,  the  fractional  change  in  C02  density  is  related 
to  the  fractional  change  in  N2  density  by  the  expression 


d[C02]  d[N9] 

- L.  =  F (a)  - .  (A31) 

[co2]  [n2] 


If  we  assume  F(a)  is  constant,  determined  by  some  appropriate  mean 
value  for  a  as  the  N2  density  changes,  then  we  can  integrate  Equation 
(A31)  to  obtain 

[C02]f  =  [C02].{[N2]f/[N2]i}F(a)  .  (A3 2) 

Thus,  the  change  in  [CC>2]  is 

<5[C02]  =  [C02]^  -  [C02]. 

=  -[C02]  (1  -  {[N2]  /[N2]  }F(a)  )  .  (A3 3) 

i  f  i 
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